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ABSTRACT 


This report presents interferometry theory for the Block II 
processor, including a high-level functional description and a 
discussion of data structure. The analysis covers the major 
processing steps: cross-correlation, fringe counter-rotation, 

transformation to the frequency domain, phase calibration, 
bandwidth synthesis, and extraction of the observables of 
amplitude, phase, phase rate, and delay. Also included are 
analyses for fractional bitshift correction, station clock 
error, ionosphere correction, and effective frequencies for the 
observables . 
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DEFINITIONS OF SYMBOLS 


a l 


= first harmonic of quantized stopping function 


A(k,a>) 


A s 

A(co) 

b 


random Fourier amplitude at frequency u received 
from direction k for an extended source 

estimated amplitude in fringe fitting 

random Fourier amplitude at frequency to for a 
point source 

bit interval, typically 250 nsec 


Sr 

c 

c k 

c mk 

d a 

Dn 

d S 

d(t) 

dn 


= retarded baseline 
= speed of light 

= relative fringe amplitude, bin k 
= model for relative fringe amplitude, bin k 
= D N /(2 W) 

= amplitude normalization 
= source diameter 

= normalized noise added to calibration signal 
= differential solid angle 


E(x, t) 

G 2 

Gj(co) 
h(tojj , t(j) 

h‘ (a)fc,t(3) 
h"(a)k,t G ) 

h c ( u k' tG) 


electric field at x, t 

quantization or fractional bitshift filter 

average of square of bandpass amplitude across W 

bandpass amplitude at station j 

initial frequency-domain fringes 

frequency-domain fringes after fractional bit- 
shift correction 

frequency-domain fringes after phase calibration 

frequency-domain fringes after amplitude 
normalization 
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k 

A 

k a 

k e 

k m 

ft 

N c 

% 

N t 

qk 

Rb ( u , v ) 
r m 

R m (x f ,o) f ) 


r(t G , ft) 
R(T£,0)f ) 
R’ ( Tf / (df ) 


Sc 

Sp 

Sp(oj) 
S«p( Aco) 


t 

t l 


apparent wave vector 

apparent wave vector from reference point in source 
proportionality constant in ionosphere phase 

model for direction of source wave vector 

lag offset of correlator 

number of BWS channels 

normalized additive noise, station k 

number of bits in fringe-stopping sum 

normalized natural source signal, station k 

brightness transform 

"calibrated" fringe amplitude 

model fringes in BWS-delay/f ringe-f requency 
domain • 

cross-correlation function 

R'(xf,cof) corrected for FFT origin offsets 

initial fringes in BWS-delay/f ringe-f requency 
domain 

correlated flux 

total source flux 

natural source power spectrum 

frequency resolution function resulting from FFT 
over lag 

true time 

value of t G for first point 


tbj 

tbk 

T c 


= shorthand for arguments in Eq. 6.7 
= bit time, station k 

= correlation interval of 25 msec (10 5 bits) 
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tck 

tG 

ii 

tG 

Tn 

t 0 

tpk 

T q 

T S 

t s 

t z 

u 

Uft(t G ) 

V 

v f 

Vk 

W 

-»• 

X 


= 4>-cal time, station k 
= common Greenwich time 
» t G ~ t]. 

= maximum integration time 
= additive noise temperature 
= true time "origin" 

= true time that corresponds to bit time tbk, 
station k 

= natural-source noise temperature 

= total or system noise temperature 

= correlator reference time 

= reference time in fringe-fitting 

= component of sky-projected baseline vector in 

direction of right ascension in units of wavelength 

= complex stopped fringes for lag fi. at time t G 

= component of sky-projected baseline vector in 
direction of declination in units of wavelength 

= fringe visibility = correlated flux/total flux 

= baseband voltage signal, station k 

= nominal recording bandwidth = sample rate/2 = 2 MHz 
= position vector in geocentric coordinates 


acfc(t) = location of station k as function of true time 

%ik = model for position of station k 

yjj = doppler-shifted frequency, station k 

Zfc = normalized analog signal at station k 

a = right ascension 

ajK = composite effect of instrumental phase shifts, 

station k 

13 = structure coordinate in direction of right ascension 
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Y 

Af s 

A k 

AN e 

Ax c 

Ax cs 

Axg 

Axr 

Axr 

Ax t 

a*4> 

A4>f 

A 4>h 

A<t>|(w n ) 

A «C 

5 

e A 

e k 

Hk 

©k 

e q 

K k 

X 

T ak 

T bk 


= structure coordinate in direction of declination 
= spanned bandwidth 

= fixed integer-bit delay introduced by correlator 

= differential electron content through ionosphere, 
station j minus station i 

= residual clock offset 

= correction required to remove effects of clock errors 
= residual geometric delay 
= residual delay in calibrated fringes 
= residual delay in uncalibrated fringes 
= residual troposphere delay 

= maximum residual phase delay rate 
= stopped fringe phase 

= "stopped heterodyne phase," station k* 

= stopped phase of n th tone, station k 
= true frequency minus nominal frequency 
= declination of source 

= fractional decrease in fringe amplitude 
= tweak frequency to remove effect of lag offsetting 
= all noise added to recorded signal, station k 
= baseband signal phase at station k 
= phase of quantization filter, Fg 
= antenna sensitivity at station k 
= wavelength of recorded radio signal 
= antenna delay* before injection point, station k 
= error in bit time*, station k 


*Use without subscript denotes difference j minus i 


x 



T BWS = final BWS delay 

x c k = error in <|>-cal time*, station k 

x^k = model clock "delay," station k 

= geometric and troposphere delay for station k 
x e = group delay due to charged particles 

xf = BWS "search" parameter in FFT 

Tg = conventional geometric delay 

T< ?mk = model geometric delay removed, station k 

xjk = instrumental group delay*, station k 

x m = T mi ~ x mj 

T MA = total model delay to be restored 

T mg = model geometric delay restored 

= retarded correlator model delay, station k 

x Mk = unretarded model delay, station k 

x s = estimated residual delay from fringe fitting 

T^k = troposphere delay for station k 

= model troposphere delay, station k 

x uk = delay between clock reference point and 4>-cal 

injection point, station k 

t x = interim geometric delay 

T<f> = final phase delay 

= estimated phase-delay rate in fringe fitting 
4 >b = phase shift due to source structure 

<j> c = oscillator phase at clock reference point 

<t> e = phase shift due to charged particles 

4>f = total fringe phase before counter-rotation 


*Use without subscript denotes difference j minus i 
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♦h 

<t>hk 

♦iK 

4>s 

<t> T (a> k ) 

Yf 

H'Sk 

CO 

CO 

C01 

to c 

(Of 

u hk 

“hk 

co k 

co k 


= total heterodyne phase shift 
= heterodyne phase shift, station k 
= instrumental phase shift*, station k 
= model phase for counter-rotation 
= phase of n th <j>-cal tone at injection point 
= estimated phase in fringe fitting 
= calibration phase at frequency co k 
= corrected fringe phase 

= residual fringe phase after counter-rotation 
= stopping phase for n th tone, station k 
= RF frequency in geocentric coordinates 

= estimated "center" frequency of 2-MHz passband 
used in fringe-counter-rotation 

= RF frequency of first frequency bin 

= nominal oscillator frequency 

= fringe frequency "search" parameter in FFT 

= true total heterodyne frequency, station k 

= estimated value for cuh k 

= RF frequency for k th bin 

= co k - c%i = baseband value for frequency co k 


“k “ «k ” W 1 

cd£ = frequency of nth <t>-cal tone 

cot = true oscillator frequency 

co z = reference frequency in fringe fitting 

co T = effective frequency for BWS delay 


co* = effective frequency for delay-rate observable 

co T <j> = effective frequency for phase-delay observable 

co<j> = effective frequency for phase observable 
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SECTION I 


INTRODUCTION 


Over the last few years, a very long baseline interferometry 
(VLBI) processor, the BLOCK II (BLKII) processor, has been under 
development at JPL and the California Institute of Technology. 

The major subsystems and their principal contributors are 
correlator software (D.H. Rogstad and H.W. Hammond), correlator 
hardware (J.C. Peterson, B. Rayher, J.W. Dillon, B. Espinoza, 
and E. H. Sigman), Tensor (M.S. Ewing and J.L. Vavrus), and 
post-Tensor software (T.G. Lockhart). 

Reading each tape channel at a rate as high as 4 Mbits/sec, the 
current version of the BLKII processor can simultaneously reduce 
tapes recorded at four stations, including cross-correlation of 
28 tape channels for all six station pairs. A single pass 
through the BLKII processor can reduce many gigabits of data to 
only four output observables: bandwidth-synthesis (BWS) delay, 
phase delay, phase-delay rate, and fringe amplitude, with each 
corrected for instrumental and processing effects. 

This report outlines the functional design of the BLKII processor 
and augments a previous analysis [1] of the BLKO system to obtain 
a mathematical description of data reduction in the new proces- 
sor. Even though the theory of data reduction in the BLKII pro- 
cessor is very similar to the BLKO analysis through the step of 
phase calibration, this report will reproduce much of the an- 
alysis and discussion in that report for completeness. As in the 
BLKO analysis, only nonrelativistic theory will be considered. 

The analysis includes all signal processing steps, beginning 
with the free-space form of the incoming radio wave and tracing 
its evolution through the station instrumentation, cross-correla- 
tion, fringe counter-rotation, transformation from the lag domain 
to the frequency (RF) domain, phase calibration, fractional— bit- 
shift correction, coherent integration over time and frequency, 
observable estimation and model restoration. In addition to 
system improvements, some of the new features in the analysis 
include: a new fringe model for least-square estimation of the 

four output observables (Section XIV); a demonstration of the 
approximate equivalence of this fringe-fitting procedure to an 
implicit phase-fitting procedure (Appendix C); formulas for 
calculating an optimal reference time and reference frequency 
(Appendix C); a method for calculating effective frequencies to 
be used in dual-band calibration of charged-particle effects 
(Appendix E) ; and a method for correcting the final observables 
for absolute clock error at each station (Section XII). 
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SECTION II 


SUMMARY OF THE BLOCK II SYSTEM 


The primary features of the BLKII processor from input bits to 
output observables are shown in Figures l to 5. Figure 1 shows 
the major steps in the BLKII processor; Figure 2 defines the 
indexing convention for station pairs; Figure 3 describes in 
more detail the processing of a pair of stations through the 
step of fringe counter-rotation; Figure 4 illustrates with bit- 
stream manipulations the steps shown in Figure 3, but for a 
single lag; and Figure 5 illustrates the structure of the data 
in the steps from cross-correlation to output observables of BWS 
delay and delay rate. 

Data acquisition for the BLKII processor is handled primarily by 
the Mark III (MKIII) data acquisition system [2], which can 
simultaneously and continuously record as many as 28 channels on 

1- inch instrumentation tape. In its most versatile form, the 
MKIII system can record each channel at any one of seven rates 
(0.125, 0.25, 0.5, 1, 2, 4, or 8 Mbits/sec). Present systems, 
however, are typically configured so they can record each channel 
at only one rate, 4 Mbits/sec. The Mark III system down-converts 
channels in upper-sideband/lower-sideband pairs, which leaves 

14 channel-pairs that can be arbitrarily divided between bands 
(e.g., S band and X band). Although the Mark III system is the 
primary system for the BLKII processor, Mark II tapes can also 
be processed. 

Even though the BLKII processor can accommodate a number of the 
potential record rates of the MKIII system, this report will 
assume that the nominal channel bandwidth is 2 MHz and the record 
rate is 4 Mbits/sec. Results for other record rates can be 
obtained by scaling according to record rate wherever necessary. 
Reading each channel at 4 Mbits/sec, the BLKII processor can 
simultaneously process tapes from as many as four stations, which 
are combined in pairs to form as many as six baselines. As 
illustrated in Figures 3 and 4, the bitstream from each station 
is offset by an a priori delay computed as the difference in 
time of arrival of the incoming radio signal at the station and 
at the center of mass of the earth. For geophysical applica- 
tions, eight lag offsets relative to the a priori delay model 
are applied to one of the two bitstreams on each baseline. (The 
spectral line option is considered below.) After multiplication 
of the multiply-offset bitstreams of a given station pair, the 
product for each lag is counter-rotated in quadrature to remove 
rapid phase variations and the result is summed over 25 msec 
(0.1 Mbits). The resulting complex fringes for each channel are 
then transformed by a discrete Fourier transform (DFT) over lag 
to the (RF) frequency domain where, for an 8-lag transform, each 

2- MHz channel is sampled at the center of each of four bins of 
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STATION TAPES 



AMPLITUDE, BWS DELAY, DELAY 
RATE, PHASE DELAY PER 
BASELINE/RF BAND/OBSERVATION 




DENOTES 6 BASELINES 
AND 2 BANDS 


Figure 1. Functional Block Diagram for the Block II Processor 
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SINGLE-STATION OPERATIONS TWO-STATION OPERATIONS 



Figure 2. Indexing Convention in the Block II Processor 
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5b. Schematic Illustration of Data Structure in the 
Block II Processor: Part II 
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0.5-MHz width. In this domain, the phase of each bin can be 
easily counter-rotated to correct for the fractional bitshift 
error caused by quantization of the delay offsets applied to the 
bitstreams. 

The fringes at this point are still corrupted by instrumental 
phase shifts and instabilities. To calibrate these effects, cal- 
ibration tones are injected during data-taking at the front of 
the instrumentation so that they will traverse the instrumenta- 
tion along with the natural source signal and will be imbedded 
in the recorded signal. The first version of a calibration sys- 
tem using tones was developed by Rogers [3]. While the bit- 
streams are being cross-correlated for each baseline, each bit- 
stream is also counter-rotated separately on the basis of tone 
model phase in order to extract the phase of each phase-calibra- 
tion tone (phase-cal tone) . The resulting tone-phase values are 
used to counter-rotate the frequency-domain fringes for each 
frequency bin in order to remove instrumental phase effects and 
thereby make the fringes coherent between RF channels. 

After these corrections, the correlator sums the fringes in each 
frequency bin over a number of time points to compress the data 
further. The integration time covered in this sum should be no 
greater than 1/12^“ the period of the highest expected residual 
fringe rate if amplitude loss is to be restricted to 1% or less. 
The separation between output fringe points will typically be 
1-10 seconds. 

After the fringes leave the correlator, the Tensor software per- 
forms a fast Fourier transform (FFT) over frequency and time, as 
illustrated in Fig. 5. First, the Tensor subsystem coherently 
sums the fringes over channels (four bins per channel) within a 
band by means of an FFT, which transforms the fringes from the 
(RF) frequency domain to the BWS-delay domain. The effect of 
this transform is to synthesize a bandwidth covering the RF band 
(e.g., 8300-8700 MHz). The BWS technique was developed by 
Rogers [4] to provide lower detection thresholds and more 
accurate delay measurements. Following the BWS FFT, an FFT over 
the time variable (e.g. , over 200 seconds) transforms the fringes 
to the phase-delay-rate domain (fringe-frequency domain). For 
each RF band, the resulting fringes are then a complex function 
of the FFT variables (i.e., the search variables of BWS delay 
and fringe frequency), a function that peaks at the actual delay 
and delay rate. (One also has the option of transforming the 
fringes over only one axis, either frequency or time, and out- 
putting the resulting array without estimating delay or delay 
rate . ) 

Signal detection consists of a search over the BWS delay and 
fringe frequency to find the maximum fringe amplitude (root-sum- 
square (RSS) of real and imaginary parts) . To improve effi- 
ciency, the range of BWS delay and fringe frequency values 
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covered in the search can be narrowed on the basis of a priori 
estimates, if desired. 

In the post-Tensor software, the final observables (amplitude, 

BWS delay, phase delay, and phase-delay rate) for each band 
(e.g., S- and X-band) are extracted by means of a multiparameter 
least-squares fit of a fringe model to the fringes observed in 
the region of maximum amplitude. The fringe model used in the 
fit is computed as a coherent sum of phasers whose time— indepen- 
dent amplitudes can be separately specified for each channel by 
the user and whose phases are a linear function of BWS delay, 
phase-delay rate, and phase delay, all referenced to a time near 
the center of the observation and to an RF frequency near the 
center of the spanned bandwidth. After estimation, the post- 
Tensor software restores the model delay removed by the corre- 
lator, thereby producing total output values for BWS delay, 
phase-delay and phase-delay rate that are theoretically free of 
processing effects. 

In spectral line applications, more than eight lag offsets are 
required to obtain adequate RF frequency resolution. For the 
BLKII processor, the number of lags per channel can be increased 
by decreasing the number of channels and using the correlators 
of the omitted channels to process more lags for the remaining 
channels. In the extreme case of one channel, 224 (28*8) lag 
offsets can be applied, which leads to a sample spacing of 
approximately 2000/112 = 18 kHz across the 2-MHz passband. 

Another BLKII feature useful in astronomy applications is the 
capability of simultaneously processing right-circularly polar- 
ized and left-circular ly polarized data. In this mode, each (RF) 
frequency channel is recorded on two tracks rather than one, so 
that a maximum of 14 channels can be recorded on MKIII tapes. 
Since cross-correlation of a given channel— pair now requires 
four combinations (R*R, R*L, L*R, L*L) rather than one, a maximum 
of seven (28/4) channels of polarized data can be cross— corre- 
lated in a single pass through the processor. With more than 
eight lags per pair, even fewer channels could be processed in 
a single pass . 


2-9 



SECTION III 


IMPROVEMENTS 


Many features of the BLKII processor represent significant 
improvements relative to its predecessor, the BLKO processor. 
These include: 

a) More (2-MHz) channels (28 rather than 8). 

b) A higher total bit rate through continuous processing 
of each (2-MHz) channel. 

c) An improved correction for bitstream-offsetting 
("retarded station") effects. 

d) An improved f ractional-bitshift correction through more 
frequent (every 25 msec) application. 

e) Better phase calibration by frequently (every 25 msec) 
applying intra-passband corrections. 

f) Reduced intermodulation and aliasing in counter-rotated 
calibration-tone phase (cal-tone phase) by quantizing 
the tone-stopping sinusoids to 128 amplitude levels 
rather than 3 . 

g) Simpler and faster delay estimation, as well as a lower 
fringe detection threshold, by coherently summing the 
fringes from intrapassband bins across all channels in 
a band. 

h) Improved capability for spectral line applications by 
cross-correlating as many as 224 lags. 

The BLKII processor could be improved in a number of ways. 
Possible improvements include: 

a) Increasing the number of stations to seven (and with 
more effort, to any number), including the associated 
21 baseline correlators. 

b) Increasing the number of phase-cal correlators from one 
to two (or four) per channel, which would allow more 
frequent (or continuous) sampling of calibration phase. 

c) Increasing the delay lines to 0.5 sec (an improvement 
required for interferometry between a ground-based 
antenna and an orbiting satellite) . 

d) Applying a digital filter to the time sequence of 
counter-rotated fringes from each frequency bin in 
order to minimize the output sample rate in time. 
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e) Improving Tensor to increase the allowed range of 

uncertainty on a priori values for baselines and source 
positions. 

Future resources will determine when these improvements will be 
implemented . 
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SECTION IV 


CLOCK MODELING 


Before proceeding to cross-correlation analysis, it will be help- 
ful to define a few time-keeping terms. First, define the fol- 
lowing thus : 

t = true time (e.g., "perfect" universal time, 
coordinated (UTC)) 

cot = true oscillator frequency (relative to true time) 

a) c = nominal oscillator frequency (e.g. , 5 MHz) 

Throughout this analysis, true time will serve as the "founda- 
tion" of timekeeping since, in a nonrelativistic analysis, it is 
the uniform time scale applicable to both antennas. 

Two "clocks" will be defined to assist in the analysis. The time 
of the first clock, the "recorder clock," is registered by the 
bits of the tape. This bit time is given in terms of true time 
by 


t 


b 


t + x b 


(4.1) 


where x^ is the error in the recorder clock. 

The second clock, the "4>-cal clock," will be established, at 
least conceptually, near the point where the 5-MHz signal from 
the station enters the phase calibrator instrumentation, as sche- 
matically illustrated in Figure 6. The phase at this reference 
point is given in terms of true frequency by 

4> c = &> t (t - t Q ) (4.2) 


where t Q is the (unknown) true time at which 4> c = 0. As implied 
by Eq. 4.2, we will assume that the oscillator is perfectly sta- 
ble so that the only imperfections in the station clock are off- 
sets in epoch and rate. If phase 4> c is measured at the reference 
point, 4>-cal time can be defined by 


t 


c 


<t> /to 
c' c 


(4.3) 


where we assume t c is zero when 4> c is zero. This is possible 
since a particular zero crossing can be defined to be t c = 0 . If 
the nominal frequency is related to the true frequency by 

co c = co t + Aco c (4.4) 
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then 4>-cal time is related to true time by 


t = t + T 
C C 


where the error in the 4>-cal clock is given by 

Aco_ 


T C CO . 


(t - t 0> - t o 


(4.5) 


(4.6) 


Usually the difference between w c and w-t is small, less than a 
part in 10 11 . 

It is conceptually and instrumental ly convenient to define a 
clock in terms of a l-pps signal that drives a counter to reg- 
ister the number of received pulses. To establish the l-pps 
signal, a particular set of positive-going zero crossings (one 
every 5 x 10 ^ crossings for 5 MHz) must be selected. (Equi- 
valently, in Eq. 4.3, a definition of 4> c = 0 must be made.) 

This arbitrary selection for the 4>-cal clock leads to an ambi- 
guity which must be resolved in clock-synchronization measure- 
ments, as discussed in Section XIII. In this report, clock 
synchronization will be specified relative to the 4>-cal clocks. 
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SECTION V 


THE RECORDED NATURAL SOURCE SIGNAL 


This section derives an expression for the analog signal by 
tracking the incoming signal from its free-space form through 
the instrumentation to its final form at baseband. Only non- 
relativistic theory will be considered in this report. 

In typical radio interferometry measurements, two antennas simul- 
taneously receive the signal produced by a very distant, compact 
natural source of random radio emissions. In this section, it 
will be assumed that the source is a point source. Effects 
caused by extended structure are derived in Appendix A and will 
be accounted for in the next section. Due to the great distance 
of the natural source, the incoming signal can be modeled by 
plane waves. Typical natural sources used in JPL applications 
emit noise with a smooth frequency distribution. A given member 
of the random ensemble of possible noise waves can be modeled as 
a superposition of plane waves of different frequencies: 


/ oo 

A(co) e^ ut k* x ^dco (5.1) 

00 


\ghere^E is the electric field at true time t and at point x, 
k = cok/c is the apparent wave vector of the plane wave, and 
A(oj) is the (random) Fourier amplitude at frequency w. All 
quantities are measured with respect to a frame (the "geocentric 
frame") whose origin is at the center of mass of the earth and 
whose axes are defined by true equatorial coordinates of date. 
Frequency co is specified in the geocentric frame on the basis 
of true time. The adjective "apparent" refers to the aberration 
effect that will be observed in a geocentric frame. For simpli- 
city, the wave is assumed to be linearly polarized. At this 
point, it will be assumed that the wave propagates in a vacuum. 
Transmission media delays will be added in the next section. 

Relative to geocentric coordinates, the electric field at a 
given station becomes 


Ej (t ) 


E(Xj ( t ) ,t) 



i[oot-ic«x . (t) ] 
e ^ do> 


(5.2) 


where Xj(t) is the location of station j as a function of true 
time t. We will assume that each antenna is a parabolic dish 
with two rotation axes, and that station location is defined to 
be the intersection of those axes. For some antennas, the axes 


5-1 


of rotation do not intersect and a correction [5] must be applied 
to the observed delay to establish an effective intersection of 
axes. This report will ignore possible axis offsets under the 
assumption that corrections will be applied in subsequent pro- 
cessing, if necessary. 

After injection, the signal passes through several filters and 
is heterodyned to baseband under the combined effect of various 
mixing signals. The phase effects of all intervening components 
can be grouped into three catagories. First, the overall effect 
of the mixing signals can be described by one total mixing phase, 
u>hjt + 4>h j ' where co^j is the total mixing frequency and rep- 

resents both a constant phase offset and all variations inmixing 
phase that are nonlinear in time. Second, all group delays after 
the injection point, including effective group delays through 
filters, can be represented by one total delay, xj j . Third, all 
remaining instrumental phase shifts can be represented by one 
total shift, <t>i j • As implied, this frequency-dependent term 
contains all phase shifts left after accounting for the hetero- 
dyne phase and the effective group delay term xjj . 

Another "instrumental" delay, x a j , will be included below to 
establish the definition of antenna location. This delay, the 
"antenna" delay, is equal to the difference in path length for 
the radio wave in its actual propagation to the tone injection 
point and its theoretical propagation to the intersection of 
axes. Since the antenna continually rotates to point at the 
source, the injection point and the intersection of axes will 
have the same geocentric velocity and there will be no differ- 
ence in the doppler shift at those two points. Thus, to include 
the antenna delay, one merely offsets the plane wave phase by 
a constant phase shift equal to the antenna delay expressed in 
terms of RF wavelength. It will be assumed in this report that 
the antenna delay is constant for all times and all sources and 
that dispersive phase shifts above the phase-calibration injec- 
tion point are negligible. If either of these assumptions is 
not valid, then additional calibrations will be required to 
remove these effects. 

Figure 6 schematically illustrates the modifications applied to 
the signal phase as the signal progresses through the station 
instrumentation from injection to the recorder. The block dia- 
gram has been grossly simplified to show only one element of 
each type, with token phase shifts and group delays. Never- 
theless, the following analysis is quite general if the var- 
iables discussed above are properly defined to represent lumped 
effects of all elements in the actual station instrumentation. 

For more details concerning the MKIII data acquisition system, 
see [ 2 ] . 
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When lumped instrumental effects are included, the doppler- 
shifted analog signal at baseband for a given BWS channel at 
station k at bit time t£,k can be written in the form: 


v„(t K1 ) 


= /“ 

Jo 


i0,. 


'k v v 'bk / ~ / A(co)G k (y k ) e d<4) + c,c ' + n k 


(5.3) 


where 


y. = w(l - k • x./c) 


(5.4) 


e k = “ (t pk ' T ak> - S - “hktpk - ‘"hk - <t ’lk ( >'k ) <5 - 5> 


In this expression, Gk is the amplitude of the system bandpass, 
nk is the additive noise, y^ is the doppler-shifted (station) 
frequency corresponding to (geocentric) frequency co, and c.c. 
denotes the complex conjugate. The time argument, tpk, is the 
true time at which the signal recorded at bit time t^k reached 
the phase-cal injection point. This time tag is computed from 
bit time by means of the relation: 


t pk " t bk T bk T Ik 


(5.6) 


where x^k is the error in bit time and xjk is the total instru- 
mental group delay from the injection point to the recorder. 

This time tag convention is established so that station locations 
can be specified in terms of true time when required in very 
accurate applications. Note that the antenna delay, x a k/ does 

not affect the time tag calculation for Xk(t) since, as mentioned 
above, the correction to the observed delay to account for the 
antenna delay (unlike the correction to the observed delay to 
account for the time tag effects of x^k and xjk) is essentially an 
additive constant for all times and sources. Thus, the antenna 
delay is an important consideration in clock synchronization but 
not in computation of station position. Because of this feature 
of the antenna-delay correction, the functions Xk(tpk) can be 
viewed as the location of the intersection of axes at tpk, the 
true time corresponding to bit time t^j • 

The baseband signal is sampled at a rate equal to approximately 
twice the baseband bandwidth and recorded in a bilevel mode. 

That is, one bit is recorded for each sample point and that bit 
is 1 if the signal is positive and 0 if the signal is negative. 
Subsequent analysis will ignore that code and act as though +1 
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or -1 were directly recorded on the tapes. A tilde over the 
voltage symbol will be used to denote two-level quantization. 

For purposes of analysis, it will be assumed that the record 
rate is 4 Mbits/sec. 

Before proceeding to cross-correlation analysis, we will reform- 
ulate the expression for the analog signal so all components are 
normalized in measurable units of noise temperature. For station 
k, let the normalized analog signal corresponding to be given 
by 


z k (t) = [' / V 3 k <t) + A* n k (t > 

where z, q, and n are the total signal, natural source signal 
(the first two terms in Eq. 5.3) and additive noise (the last 
term in Eq. 5.3), respectively. We will assume these quanti- 
ties are all normalized so that 




(5.7) 


2 2 2 
<z > = <q > = <n > = 1 


(5.8) 


where <> denotes an ensemble average over all possible noise 
forms. The quantities T S , Tg, and T n are the (average) noise 
temperatures for the total signal, the natural source, and the 
additive noise, respectively. For the continuum sources typi- 
cally observed in geophysical applications, it will be assumed 
that the source temperature is fairly constant across the (2-MHz) 
passband so that one average temperature can describe the whole 
passband. (The noise power per passband is actually kT s W, but 
the factor kW drops out in the ratios.) As we shall see, the 
normalized form in Eq. 5.7 is more convenient for cross-cor- 
relation analysis. The noise temperature of the source can 
be obtained from source flux through the relation 


T qk " K k S F 


(5.9) 


where Sp is the total source flux in janskies and is the 
antenna sensitivity station k in kelvins/jansky . 
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SECTION VI 


CROSS-CORRELATION 


The BLKII correlator aligns the bits on the two tapes, multiplies 
them bitwise and fringe-stops (counter-rotates) the resulting 
fringes, as illustrated in Figs. 3 and 4. This section develops 
a detailed mathematical model for the cross-correlation function, 
the expectation value (ensemble average) of the product of the 
aligned signals. The expectation value is analyzed to remove 
statistical fluctuation in both signal and noise and thereby re- 
veal the average behavior of the cross-correlated signal. Only 
one frequency channel will be explicitly analyzed even though the 
other channels are processed in parallel. 

For station k, the unretarded delay model is computed as 


- 


m 


3 W t) 


m / . % 

- T tk (t> - 


m / . \ 
T ck (t) 


+ A 


k' 


k = i or j (6.1) 


in which the first term is the geometric delay relative to the 
earth's center of mass, expressed in terms of the modeled values 

for station location x m k and apparent wave direction k m (with 
both in geocentric coordinates). The term x™k models the tropo- 
sphere delay (« l/sin EL), while T$k models the station clock 
error (including offset and rate). The delay A^, which is speci- 
fied below, represents one of two fixed integer-bit offsets 
required to center the fringes in the range of lag offsets. 

For station k, the retarded delay x m fc(t) used to offset the 
bitstream is computed as 


T mk (t) " T mk (t - T mk> 


( 6 . 2 ) 


u 


where x ^ is obtained by iteration, beginning with x mk = x m ^. 

As will be shown in Section VII, this method for retarding the 
model delay compensates for the effect of bitstream offsetting 
on the time tag of geometric and tropospheric components of the 
actual phase. 


To obtain the total bitstream offset, the correlator rounds the 
retarded model delay in Eq. 6.2 to the nearest bit, which will 

be denoted by a tilde over x m ^. In Fig. 3, the delay x^ refers 
to the first three terms in Eq. 6.1, while the delay x m includes 

the fixed delay A^ following x^. 

The correlator offsets the signals recorded at the two stations 
and multiplies them, as illustrated in Figs. 3 and 4. The ensem- 
ble average of this product, the cross-correlation function, is 
defined by 


6-1 



r<t G .H) S <V i (t M - ? mi (t G ) - lb) V.(t bj - ? mj (t G » > (6.3) 

where Vk(tfck) is the recorded (two-level) voltage at bit time t^k 
at station k; x m i and x m j are the quantized model delays used to 
offset the two bitstreams; and % is any one of eight lag offsets 
(H = 0 to 7) applied to station i. For a record rate of 4 Mbits/ 
sec, the bit interval b is equal to (4 Mbit/sec) -1 = 250 nsec. 

The quantized model delays are obtained by rounding to the 
nearest bit the delay values computed from Eq. 6.2 at time tg. 

As mentioned above, the delays x m i, x m j each include a fixed 
delay, Ak» applied by the correlator. The fixed delay is equal 
to 2 bits for station i and 5 bits for station j . Based on these 
values for the fixed delays, the peak lag is nominally Q. = 3 . As 
explained in [1], p. 69, the correlator applies to the bitstream 
of station i a pairwise vernier delay equal to (-1, 0, or + 1) 
bit in order to compensate whenever single-station rounding is 
inappropriate for pairwise offsetting. The vernier delay is 

implicitly contained in x m i in Eq. 6.3. Note, however, that 
the fixed delay and vernier delay are applied together in actual 
correlator operation, as shown in Fig. 3. 

The argument tg is the common Greenwich time defined by the bits 
on the two tapes. That is, the correlator implicitly enforces 
the assumption that 


t 


bi 



(6.4) 


and uses tg as the time tag. Corrections are applied in 
Section XII to account for time tag errors introduced by this 
assumption. 

For a small correlation amplitude (r « 1), the van Vleck rela- 
tion [6] gives the cross-correlation function in Eq. 6.3 in terms 
of the analog signals: 


- 5 


2 < V i (t bi 


*mi - * b) VVi - T pi» 


•J<v? > < v?> 


(6.5) 


(In typical VLBI measurements, one finds r < 0.2.) With the nor- 
malized representation of the preceding section, the expectation 
of the normalized voltage product in Eq. 6.5 becomes 


I 
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( 6 . 6 ) 


< v i(t bi ) v i (t bi )> 

yj < vf > < v? > 


<Zi(t 



z j( t 



Tg- < q i (t bi ) q j (t bj ) > 


(6.7) 


since the additive noise nj is uncorrelated between stations. 
The abbreviated notation for the arguments is defined by 


*bi E *bi - ’'mi - lb 


( 6 . 8 ) 


t bj t bj ^mj 


(6.9) 


Thus the problem is reduced to computing the correlation between 
the normalized natural source components (q^) for the two sta- 
tions, with the overall amplitude scaled by a multiplicative sys- 
tem temperature factor. Note that this factor times 2/tt gives 
the amplitude of the cross-correlation function for perfectly 
correlated natural source components (i.e., <qiqj> = 1.0). 

We will now present on the basis of a plausibility discussion 
the final form for the cross-correlation function. A complete 
derivation is given in Appendix A, including the effects of 
source structure and the random nature of the source signal. 

The natural source signal, q, is given by the first two terms in 
Eq. 5.3 but is normalized to have unit power. Thus the product 
q^ qj can be represented by a double integral over frequency 
where the integrand includes the sum of two products of the form 
A((o)A*(w') and A(co)A(a) ' ) . After the ensemble average, only the 
product <A(w)A*(w’)> has a non-zero value and this occurs as a 
delta function at id = io'. (It will be assumed that the signal 
from the natural source is stationary so that the amplitude at 
one frequency is uncorrelated with the amplitude at another fre- 
quency. ) Thus, the double integral over frequency collapses to 
a single integral over frequency with one term in the integrand. 
If we let Sp(oj) denote the source power spectrum arising from 
<A(co)A*(o))>T the cross-correlation function becomes 


r(t G ,a) 


2 ^m 1_ 

v 2 D N 
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00 


Sp(co) G i (y i ) 
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i<J> f 

e dco + c . c . 


( 6 . 10 ) 
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where 


r = v c 
m f 


T . T . 

- 31. -31 

T . T . 
si s] 


( 6 . 11 ) 


- 


= f S_(w) G?(y. )da) f S (co) G?(y.)dw (6.12) 

^ n ^ JJ 


*f ' ♦h - “hj *mj + “hi (? mi + 8b) 


+ k a ' [*j (t pj - V - x i <t: pi - T mi - * b) 
+ U [x b + Tj + x a - x mi + x mj - lb] 


+ *j + 4>„ + 4> e 


(6.13) 


in which 


♦h = '“hjhpj + ♦hj’ - ‘“hi^i + *hi> 


(6.14) 


and x^, xj, x a , and 4>j are differences between the corresponding 
station quantities in the order j minus i (e.g., xj = xjj - xii). 

The instrumental phase shift, 4>j, is the difference of the 
instrumental shifts at the two stations, with each given an 
argument offset to account for doppler shifting: 


+i = 'hj ( yj ) - 4, n ( 3'i ) 


(6.15) 


As noted in Section V, represents nonlinear terms in the 
phase- frequency response of the instrumentation, while the delay 
xj represents linear terms. The linear terms will be accounted 
for when the "time tag" correction is applied (see Section XII). 
The current phase-calibration procedure will remove phase effects 
to first approximation but does not account for doppler shifting 
in the nonlinear terms. If all channels have matched phase- 
frequency responses, this doppler effect would essentially cancel 
in the BWS delay. This potential error source will not be 
treated in this report, but it might be large enough in some 
applications to require correction. 
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The amplitude terms in the cross-correlation function in Eq. 6.9 
can be explained as follows. The factor of 2/ir is the amplitude 
loss due to two-level sampling (see Eq. 6.5). The factor of 1/2 
is a consequence of complex representation of the fringes and 
gives the proper definition for the power spectrum Sp (i.e., by 
cancelling the factor of 2 gained when the complex conjugate in 
Eq. 6.10 is added to obtain a real number for the cross-correla- 
tion function) . Normalization of the signal in terms of system 
temperature generates the factor r m , as suggested by the deriva- 
tion of Eq. 6.7. Included in r m is the fringe visibility Vf, 
which accounts for loss in fringe amplitude due to the resolving 
power of the interferometer. Fringe visibility can be expressed 
as 


v f = S c /S p (6.16) 


where S c is the "correlated" flux of the source and Sp is the 
total flux. The fringe visibility and a brightness transform 
phase shift $3 are two of the three modifications to the final 
cross-correlation function when the analysis is generalized to 
include extended sources as shown in Appendix A. The third 

modification is to change the point-source wave direction k to 

k a , a wave direction defined relative to a somewhat arbitrary 
reference point within the extended source. 

2 

In the denominator, the quantity Djj is the product of the total 
source power from the two stations and arises from the normaliza- 
tion of the natural source components at each station (that is, 
from the requirement <q?> = <q?> = 1). Note that, as one would 
expect, the total source power at each station is equal to the 
integrated power across the passband where the power at each 
frequency is given by the spectral power of the natural source 
multiplied by the square of the bandpass amplitude. 

Within the integral in Eq. 6.10, the various terms can be 
explained as follows. Since a component of the natural source 
signal at one frequency is uncorrelated with a component at any 
other frequency, the correlated signal can be obtained by multi- 
plying each frequency component at one station by the correspond- 
ing component at the other station and then summing this product 
over all frequencies. By corresponding frequency components we 
mean the components with the same frequency in geocentric coordi- 
nates. Due to the earth's rotation, a given geocentric frequency 
will be doppler-shifted when measured by station instrumentation. 
Those shifts are accounted for by the doppler-shifted arguments 
yjj in the bandpass functions and in the phase shifts 4 >ik- 
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Since interferometers are differencing devices, the phase 4>f in 
Eq. 6.13 is the difference in the f requency-specif ic phase for 
the two stations, each with the offset in time introduced by 
correlator delays. All of the phase terms in Eq. 6.13 can be 
derived in this manner except for three which involve special 
considerations. First, as mentioned above, the brightness trans- 
form phase 4 >b is a consequence of source structure as shown in 
Appendix A. Second, the quantities x^ and xj arise from the dif- 
ference tpi - tpj , which, through the use of Eqs . 5.6 and 6.4, 
can be shown to equal 


t 


Pi 



T b + X I 


(6.17) 


Third, a phase term 4> e has been added in Eq. 6.13 to account for 
differential ionosphere effects, with the station difference in 
the order j minus i . 

Introduction of an "interim" delay, x x , simplifies the expres- 
sions for the fringe phase in Eq. 6.13: 


*f m +h - “hj T mj + + lb) 

+ .“ [ T x + T b + T I + T a - *mi + Vj ~ 11 
+ <>1 + ts + *e 

The interim delay is defined by 

T x E T di (t pi - f mi - ‘ b) - T dj (t pj - Vj> 


(6.18) 


(6.19) 


where 


k • x.(t) 

T dk (t) " c T tk (t) ' k = i or j (6.20) 


Note that the definition of interim delay has been expanded here 
to include the troposphere delay x^, which had been neglected 
to this point. 

The phase rate 4>f can reach about 20 kHz at X-band. Therefore, 
before the data can be compressed by summing over bits, the 
fringe phase must be counter-rotated to a very low rate 
(< 0.1 Hz). The next section discusses that procedure. 
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SECTION VII 


FRINGE COUNTER-ROTATION 


In most observations of natural sources, the system noise is much 
greater than the correlated amplitude so that the signal-to-noise 
ratio (SNR) for a single bit product (Eq. 6.3) is very small 
(0.0002 - 0.2). Thus many bits of data (~10 3 - 10 9 ) must be 
collected to raise the SNR to a usable level on the integrated 
signal. However, due to the high phase rate that the fast 
fringes typically possess (as high as 20 kHz at X-band) , the data 
volume cannot be compressed by directly summing the fringes over 
time (bits). This section explains fringe stopping, a process 
that counter-rotates the fringes (i.e., digitally heterodynes 
the fringes to a low frequency) and thereby allows successive 
sums over time with each sum compressing a correlation interval 
of 25 msec (10 5 bits) into one point (stopped fringes or 
correlation coefficient). In addition to permitting a massive 
data compression, counter-rotation can, with modification, 
partially remove the effects of model-delay quantization and of 
bitstream offsetting. 

To avoid significant losses in output fringe amplitude, the model 
phase used to counter-rotate must be constructed with sufficient 
accuracy to insure that only a small fraction of a cycle is tra- 
versed by the stopped fringes over the total integration inter- 
val. (As indicated in the next two sections, the stopped 
fringes are ultimately summed over a time interval that might be 
of the order of 1-10 sec.) For example, if the geometric delay, 
the component with the greatest uncertainty, is known to a part 
in 10 6 (e.g., a 0.2 arcsec position error), then the stopped 
phase will change by 20 kHz»lO -6 «l sec =0.02 cycle across the 
integration interval when the fringe rate is at its maximum 
value of 20 kHz at X-band and the integration time is 1 sec. 
Amplitude loss due to a 0.02 cycle change in phase is only 0.07 
percent. It will be assumed that amplitude loss due to the 
inaccuracy of the fringe model is negligible. 

The model phase used for fringe counter-rotation is identical in 
form to the BLK0 model and is given by 


W - <“hj 


+ 0 ) 


- %i> (t G * V - *hj Vj ( V + *hi WV 

[WV ' WV ‘ WV + T ml <t G>] 


(7.1) 


where oihk is the best estimate for the total mixing frequency, 

u hk' anc * ® is close to the centroid frequency of the 2-MHz pass- 
band. The time t s is a reference time adopted at the correlator 
and usually marks the start of the observation. Note that the 
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delay for station i is set equal to the "2. = 0" value, which 
means the same counter-rotation phase is used for all lag offsets. 
The variation of fringe phase with lag (i.e., in x x ) will be 
removed below. 

To counter-rotate, the voltage product in Eq. 6.3 is multiplied 
at each bit by a phaser based on Eq. 7.1 and summed over the time 
points (bits) within a correlation interval, producing one com- 
plex value for that interval: 


U d (t G ) T mi ab) ^j (t G T mj * e 


-^m (t G> 


(7.2) 


b G 


This process is illustrated in Fig. 3. The counter-rotation sum 
covers an interval of 100,000 bits (25 msec), during which the 

bitstream offsets x m i, x m j remain constant. The time argument 
tQ now refers to the time at the middle of the sum interval. 

The BLKII correlator, like the BLKO correlator, approximates the 
fringe counter-rotation sinusoids with three-level (-1, 0, +1) 
functions (see p. 19 in [1]). We will assume that the only mod- 
ification introduced in the stopped cross-correlation function by 
this quantization is a small increase in amplitude, which will be 
applied below. Fringe SNR is reduced about 4 percent by this 
quantization. A simplified example of counter-rotation in terms 
of bits is shown in Figure 4 for an artificially short integra- 
tion interval of 8 bits. 

By combining Eq. 7.2 with the expectation value in Eq. 6.3, one 
finds that the expectation value for the complex stopped fringes 
is given by 


< VV> 



r(t^, a) 


~ i ^ro (t G ) 


(7.3) 


When Eq. 6.10 is substituted for r(t, &), two terms result: a 

sum-frequency term and a difference-frequency term. For suffi- 
ciently high fringe rates ($f >> 100 Hz), the sum-frequency term 
essentially averages to zero in the sum over time, leaving only 
the difference term: 
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< u l <t G ) > = 5 ^ Lj I S p <U) °i (y i > G j <y j > ® d “ (7 ’ 4) 

t G 0 


where the stopped phase is given theoretically by 


A +f = <t> f - <t> m 


(7.5) 


A4> h + wAt r + 4> t + 4> b + 4> e 


- (ca - oo^) fi-b + (to - to) (x m - x m ) (7.6) 


for which we have made the definitions 


lit 

H 

< 

x x + x b + X I + T a- T m 

(7.7) 

3 

III 

T ini <t G ) ' T mj (t G ) 

(7.8) 

3 

III 

x mi* t G* “ x mj* t G^ 

(7.9) 

III 

XX 

<J 

4> h - (i hj - “ hi > <t G - t s ) 

(7.10) 


Two terms of the form (o>hj - co^ j ) x m j that are of the order of 
10 -4 cycle or less have been omitted. 

The delay offsets (*f m i and ?mj ) applied to the two bitstreams 
introduce implicit effects in fringe phase that enter primarily 
in the time dependence of the actual geometric and tropospheric 
delays. One can readily demonstrate that model delay and phase 
have been constructed to approximately remove these effects. 
Expand both the interim delay in Eg. 6.19 and the model delay in 
Eqs. 6.2 and 7.8 about their time arguments and insert these 
expansions in the residual delay in Eg. 7.7 to obtain 
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4t R = T di (t pi ) ' <T mi + lb)x di • T dj <t pj ) + Vi T dj 

“ x mi (t G> + x mi x mi + x mj ( V “ V) x mj + x b + T I + x a 


(7.11) 


By design, all but one of the first-order terras in Eq. 7.11 
cancel if one neglects the small difference (~10 -12 ) between 

^dk an< * ^mk an< ^ n e < ?i e cts the quantization error (i.e., 

x d k * ( T m k “ ^mk^ ~ 1,6 * 10_6 * 125 nsec = °- 2 psec). With 

these approximations, one obtains a residual delay independent 
of correlator offsets except for a small term dependent on 2.: 


At R = 


T di (t P i> 


x £i ( V - T dj <t pj> + x mj (t G> 


- T di tb + T b + T I + T a 


(7.12) 


Removal of the 8.-term will be discussed in the next section. 

The remaining differences in arguments due to clock synch errors 
(i.e., tQ vs. tpi or tp j ) cause very small delay errors that are 
considered in Section XI I. 

The theoretical expression for the stopped fringes can be re- 
written by interchanging the sum and integral in Eq. 7.4 to 
obtain 


< U|l <t G » - N t ^ 

N 





s P <“> 




it Y f -(co-a) hi +to? di )Sb] 
e dtu 


(7.13) 


where the "quantization filter" ( "f ractional-bitshift filter") is 
defined by 
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Fg(w-oj) = 


i( CO-03 ) (t - T ) 

m m 


N 


(7.14) 




and where the residual fringe phase is given by 


i|/ f (co,t G ) = A4> h + + 4 > b + 4> e 


(7.15) 


in which the residual delay is defined by 


4t R - T di (t pi > ' T ni <t G ) - T dj (t pj > + T mj 


+ T, + T t + t, 

b I a 


(7.16) 


N t is equal to the number (10 5 ) of time points (bits) in the sum. 

The factor ai = 1.176 has been inserted to account for the 
amplitude increase caused by quantization of the stopping sinu- 
soid, as discussed above. The sum over time has been applied 
only to the " quant i zed-del ay phaser" from the last term in 
Eq. 7.6, since the time variation of that factor will dominate 
when an accurate # model delay is used. That is, the rate of the 
residual phase ($f) can be made very small (< 0.02 Hz) so that 
the stopped phase i|/f will change by only a small amount 
(< 1 mcycle) over a sum interval (25 msec). On the other hand, 

A A 

the phase (co - u))(x m - can change by as much as (to - w) 

* m At < 1 MHz* 2 ysec/sec* 25 msec = 50 mcycle over a sum 
interval for a 2-MHz bandpass. As explained in Appendix D of 
( 1 ] , the quantization filter Fq can easily be calculated on the 
basis of delay quantization applied at the correlator. Excep- 
tions to the above assumption of small residual rate can occur in 
"initialization" runs involving a poorly known baseline or source 
position. However, this will occur only for the first pass 
through the correlator since the results of that pass can be used 
to improve the a priori delay for all subsequent passes involv- 
ing that baseline or source. 
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SECTION VIII 


TRANSFORMING TO THE FREQUENCY DOMAIN 


The lag-domain fringes in Eq. 7.13 are corrupted by instrumental 
phase shifts and by the fractional bitshift error resulting from 
quantization of the bitstream offsets. These effects can be 
easily removed through counter-rotation if the fringes are first 
transformed from the lag domain to frequency domain. Based on 
the assumption that the a priori delay model is accurate enough 
to place the maximum fringe amplitude near the center of the lag 
range, the fringes are transformed by means of a discrete 
Fourier transform (DFT) as follows: 


VI 


MW 


-E 

a=o 


W 


K to k +e k) Q 'b 


( 8 . 1 ) 


where Ng_ (= 8) is the number of lags and where co£ is any one of 
four baseband frequencies and is a tweak (doppler) fre- 
quency inserted to remove the effect of lag offsetting in the 
argument of the interim delay. The four baseband frequencies, co£, 
are placed within the 2-MHz passband approximately as follows: 
0.25, 0.75, 1.25, and 1.75 MHz , where the sign is positive (nega- 
tive) for the upper (lower) sideband. To compute the corre- 
sponding RF bin frequencies ((%), one must add the total mixing 
frequency (<%i) at station i to these baseband frequencies. 

Details concerning the implementation of the DFT in Eq. 8.1 will 
not be presented in this report. That implementation involves a 
preliminary frequency-shift operation to properly place the 
final bin frequencies and to apply the correction, followed 
by a standard FFT over lag. In subsequent Tensor processing, it 
is assumed that all bin frequencies across the RF band fall on 
a uniform grid with 0.5 MHz between points. Thus, the L.O. fre- 
quency for each channel and the frequency shift mentioned above 
must be properly chosen to produce such placement of bin frequen- 
cies. If phase closure is desired and if the L.O. frequencies 
are slightly offset from one another around the baseline triangle, 
the frequency shift operation can be adjusted to produce the same 
bin frequencies, in RF terms, for all baselines. 

The tweak frequency is given a value that will approximately 
remove the phase term in Eq. 7.13: 


e k = <“k + 


“hi> 


mi 


( 8 . 2 ) 
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in which x m i is the model delay rate for station i at the center 
of the 25-msec interval. In practice, one can replace + 
oihi in this expression with the centroid frequency, co, with 
negligible loss of accuracy. For analysis purposes, however, 
the more exact expression in Eq. 8.2 will be used. The effect 
of this correction is to shift the output frequency spectrum by 
a small amount (< 20 kHz at X-band for earth-fixed antennas). 

The expectation value of the DFT is given by 


V 1 


i(co'+e. )Slb 

<h(o)k . tg)> = / X u Jt^6^ e 


1*0 


(8.3) 


where the sign convention for frequency is as described above. 

By substituting the expression for <u^> found in Eq. 7.13 and 
interchanging the sum and integral, one obtains 


0 ° 

■<h«*,t G )>- N t f s p (») G i (y i > G^y.) F q S T (A M ) 

W r\ 


where the convolution function is given by 


iy f 

e doo 


(8.4) 


S t (Aco) = e 


iAu(N^-l )b/2 sin (AwN a b/2) 
sin (Aa)b/2) 


(8.5) 


and the difference frequency by 


Aw E [co]j + (wic + uhiJxmi] - [u 1 + (cu' + cohi)Tdi] 


( 8 . 6 ) 


in which 


w ' 


w 


u 


hi 


(8.7) 
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The convolution function is centered on each frequency bin and 
represents a filter with a null-to-null width of 2/(N^b) -1 , 
which, for an 8-lag transform and a 4-Mbits/sec sample rate, 
equals 1 MHz. Because of this convolution, the output phase and 
amplitude (i.e., h(co£)) for a given DFT frequency co£ will not be 
precisely equal to the input cross-power spectrum at that fre- 
quency but will be equal to a weighted average of that spectrum 
around that frequency. As analyzed in Appendix E of [1], the 
phase discrepancies due to this effect are surprisingly small 
for an 8-lag transform. Given an "ideal" MRKIII bandpass shape 
(bandpass B in Table El of [1]), the difference between the phase 
of the convolved spectrum and the cross-power phase is less than 
4 mcycles for each bin and less than 0.1 mcycle for the composite 
phase for the whole bandpass. Discrepancies in amplitude, on the 
other hand, are as high as 6% for individual bin amplitudes and 
1.5% for the composite amplitude, where both percentages are 
relative to the maximum bandpass amplitude. 

Since the phase errors introduced by the limited-lag transforms 
are generally not significant for the present system, let the 
number of lags approach infinity for the purpose of approximately 
determining the form of the transformed fringes. As the number 
of lags approaches infinity, the convolution function becomes 


Lim S T (Ato) = 4irW S(Aco) 

N a -*x> (8.8) 


where 6(x) is the delta or impulse function and W is the record- 
ing bandwidth (2 MHz). As indicated by Eq. 8.6 for Aco = 0, the 
delta function enforces the following relation: 


• A • 

0) 1 + ( CU ‘ + (0, . )x,. = CO.' + ( CO.' + CO. .)x • 

v hi' di k ' k hi' mi 


(8.9) 


Since the lag-offset effect is small (cox^i < 20 kHz at X-band) , 

-12 

since co^^ is nearly equal to co^ (to 10 ), and since the delay 

model is accurate, we readily see that the delta function chooses 
the following RF frequencies: 



+ co. 


hi 


( 8 . 10 ) 


That is, the DFT outputs frequency-domain fringes at four fre- 
quency points whose RF values are equal to the four baseband 
values given in the first paragraph of this section plus the 
total heterodyne frequency of station i, co^i • For an infinite 
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number of lags, the delta function in Eq. 8.8 can be substituted 
for St in Eq. 8.4 to obtain the output frequency-domain fringes: 


<h( M . t )>- N * a r gBi^ - gL (Y .l l ) .. ?J (Y 3fc ) F e i »f ( “k* t G ) 
' n ' M k 1 N t ir a l r m D« *q e 


( 8 . 11 ) 


where 


D N 

D = 


A ~ 2irW 


( 8 . 12 ) 


Yik = - k • Xi/c) 


and where 

ff = A< *“h + "k AT R + +i + +b + ♦e 

in which Atr is given by Eq. 7.16, the quantization filter 
Fg by Eq. 7.14, Djj by Eq. 6.12, and r m by Eq. 6.11. 

In the following analysis, we neglect the small discrepancies 
resulting from the limited-tag transform and use Eq. 8.11 to 
describe the fringes in the frequency domain. If necessary, 
the transformation discrepancies could be corrected to a large 
extent through modeling, but such measures will not be treated 
in this report. 

When phase closure observables are formed in source structure 
measurements, processing effects must "close" (sum to zero) 
around the baseline triangle. Closure of processing effects 
requires antisymmetry in residual phase with respect to inter- 
change of stations. Note that residual phase in Eq. 8.14 is 
antisymmetric with respect to interchange of the i and j sub- 
scripts except for one difference implicitly contained in 
Atr, namely Aj - Aj, which is always equal to +3 lags. To enforce 
antisymmetry in this difference, the fringes can be 
counter-rotated at each bin frequency by a phaser of the form 
-3iw^b 

e , which, in effect, changes the value of Aj from 

5 to 2 bits. It will be assumed in subsequent analysis that 
such a rotation has been performed and that the value of Aj 
is changed to 2 bits in Eq. 6.1. 

The DFT over lag has cast the fringes in a form in which the 
phase effects due to instrumentation and f ractional-bitshift 
error are exposed as phase shift in each frequency-bin phaser. 

In this form, these unwanted effects can be removed through 
counter-rotation, as explained in the following sections. 


(8.13) 


(8.14) 
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SECTION IX 


FRACTIONAL BITSHIFT CORRECTION 


When the bitstreams are offset by the model delays, the offsets 
are rounded to the nearest integer bit. This rounding intro- 
duces a fractional bitshift error that ultimately shows up in 
Eg. 8.11 as an unwanted phase shift residing in the quantization 
filter, Fq (see Eq. 7.14). This section outlines the removal 
of the fractional-bitshift error in fringe phase. 

Since the BLKII processor does not change the quantized delays 
(bitstream offsets) during a 25-msec counter-rotation sum, the 
quantization filter in Eq. 7.14 becomes 


F g<“k - “> * 


sin[ (u)£ - t»>? m T c /2] ie 
e 


K - “) ? 


(9.1) 


x T /2 
m c 


where Tp is the sum interval and x m is evaluated at the center 
of the interval. For an integration time of 25 msec, a 2-MHz 

bandwidth, and an intercontinental value of 2 ysec/sec for x m , 
the sinx/x amplitude loss at the edges of the 2-MHz bandpass will 
be approximated by 0.5%. This small amplitude loss will be 
ignored in subsequent equations. The "fractional bitshift" phase 
in Eq. 9.1 can be derived from Eqs. 7.8, 7.9, and 7.14: 


0 ss (to, — w) X .(t„) - T .(t„) - X .(t„) + X .(t„) 

q v k ' |_ mi G m;) v G mi ' G mj v G_ 


where co is the "center" frequency used in Eq. 7.1 and where t@ 
now refers to the center of the 25-msec interval . 

This phase correction is the conventional fractional bitshift 
error (i.e., roundoff error) at the center of the 25-msec corre- 
lation interval. Since these phase terms are accurately known, 
they can be counter-rotated out of the fringes given in Eq. 8.11. 
The expectation value of the resulting fringes is given by 


2 S D (u k ) G i (y ik } G i (y ik ) 

<h ’ (to ,t )>x N - a r p K — — 1 — ^-e 1 K ° 

kG tirlm D 


(9.3) 


where n>f, given by Eq. 8.14, is now the total phase. 
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SECTION X 


PHASE CALIBRATION 


In radio interferometry measurements, the observed phase is cor- 
rupted by unknown and unstable shifts due to instrumentation. 

Such phase effects can degrade the accuracy of geophysical mea- 
surements and complicate measurements of clock synchronization. 
Most of these instrumental effects can be removed through the 
use of a phase calibrator. The commonly used approach is to 
inject, at a point near the front of the instrumentation, a cali- 
bration signal consisting of a set of tones, equally spaced in 
frequency and derived from the system frequency standard. When 
properly applied, this technique will both correct for system 
phase instabilities and allow absolute clock synchronization. 

To extract the calibration tones, the correlator counter-rotates 
the bitstream from each station with constant-frequency phaser 
models. The current version of the BLKII processor time-multi- 
plexes between the tone models, dwelling for 25 msec on each. 

The phase of the counter-rotated tones is extracted with an 
arctan operation. The resulting tone phases contain all instru- 
mental effects between the injection point and the recorder and 
thus can be used to remove those instrumental effects. 

Two phase calibrators are currently in use: the Rogers version 
and the Sigman version. The Rogers phase calibrator [3] gener- 
ates tones that are spaced by 1 MHz. We will assume that, when 
this version is used, the total mixing frequencies have been 
chosen so that the two tones in each 2-MHz passband have loca- 
tions suitable for interpolating or extrapolating phase to other 
frequencies in the passband. Before interpolation, one should 
carry out additional calibrations to obtain an estimate for the 
phase-frequency slope of the tone phases so that the slope ambi- 
guity (n x l nsec) due to integer-cycle ambiguities in tone 
phase can be removed (see Appendix G in [l]). if absolute cali- 
bration of the tone phase slope is not available, then care must 
be taken to see that each channel and each observation is reduced 
on the basis of the same ambiguity- induced slope error. Result- 
ing BWS delay values will then all be in error by the same addi- 
tive constant. In this case, clock synchronization measurements 
will be corrupted but geophysical/astrometric measurements will 
be unaffected. 

When the Sigman phase calibrator is used, the tone spacing can 
be set to any one of a number of possible values (see Appendix 
B) . We will assume that the selected spacing with the Sigman 
calibrator is 0.5 MHz and that the total mixing frequencies are 
chosen so that each of the four DFT frequencies (see Section 
VIII) for any given 2-MHz channel is nearly equal to a tone 
frequency. In this case, interpolation is not necessary since 
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the tone phases are extracted at the desired frequencies. How- 
ever, if clock synchronization is a goal, the a priori delay, 
including clock offset, still must be known well enough to 
resolve the ambiguity (n»2 psec) in the "delay-resolution 
function" (see Section XIII). 

In the following analysis, it will be assumed that phase cali- 
brators are available at both stations and that all necessary 
steps have been taken to produce tone phase values at the 
desired passband (i.e., the DFT) frequencies. 

The calibration phase at frequency within a passband is 
calculated as the difference in the tone phase that the 
correlator obtained for two stations: 

4> T (cj k ) = A4> i T (o) k ) - M>j T (co k ) (10.1) 


Based on Eq. B.8 in Appendix B, this phase is theoretically 
equal to 

* T( V = “k (T b ■ T c + T u + T I ) + *1 + Acp h 


where x u is the station-differenced cable delay from the 4>-cal 
reference point to the injection point and t c is the station- 
differenced synchronization error of the phase-cal clocks, with 
the differences in the order j minus i. For simplicity we have 
assumed that all dispersive phase shifts in 4>ik i- n Eq. 5.5 occur 
after the tone injection points so that they appear identically 
in the tone and fringe phases. If significant dispersive phase 
shifts occur before the injection point, separate calibrations 
will have to be applied. 

When the fringes in Eq, 9.3 have been counter-rotated with the 
tone phase in Eq. 10.2, the expectation value of the fringes 
becomes 


/h „, t rj 2 S p (c0 k ) G i (y ik ) G j (y ik ) 1 ' v c (a) k' t G ) 

< h ^k' t G >~ N t * a l r m Dl e 


(10.3) 


where the corrected phase is given by 


WV = ’MvV “ + < “k > 


(10.4) 
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Based on Eqs . 7.16, 8.14, and 10.2, the corrected phase is 
theoretically given by 


f/VV = "k 4T R + VVV + +e (< VV 


(10.5) 


in which the residual delay is now given by 


At R = T dl (t pi ) " X mi <t G ) ^ W'pj’ + T mj (t G> 


+ T c - T u + T a 


( 10 . 6 ) 


The frequencies 0 % are defined by Eq. 8.10. 

The fringe phase is now theoretically free of instrumentation 
effects introduced between the phase-calibrator injection point 
and baseband. 

After these corrections, the fringes are summed over a number of 
time points to further compress the data. If the maximum allowed 
amplitude loss due to incoherency is 1%, the integration interval 
for that sum is set by the user to a value no greater than l/12 th 
the period of the highest residual fringe rate. For example, if 
the maximum residual fringe rate were 10 mHz, then the integra- 
tion time would be 100 sec/12 = 8 sec. In the following equa- 
tions, it will be assumed that the loss in amplitude resulting 
from the sum is negligible (e.g., < 1%). Negligible phase error 
is introduced by this averaging. The form of the fringes is 
still given by Eq. 10.3 except that N t now is equal to all of 
the bits summed in the integration interval and t q refers to 
the center of the interval . 
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SECTION XI 


FRINGE NORMALIZATION 


This stage in the processing provides a convenient point to 
normalize the fringes on the basis of the number (N t ) of bits 
in the sum interval. The current version of Tensor provides two 
options: a) no normalization, which makes the fringe amplitude 
proportional to N t , and b) the conventional normalization of 
dividing by N t , which makes the fringe amplitude independent of 
N-£. The following formulation will reflect both normalization 
approaches . 

After normalization, the fringes in Eq. 10.3 can be cast in the 
form 


< h c (u) k ,1: G^> 


S c c k 


^c^k'V 


(ll.l) 


where the phase is given by Eq. 10.5. Eqs. 5.9, 6.11 and 6.16 
have been used to rewrite the amplitude in Eq. 10.3 in terms of 
correlated flux S c and relative bin amplitude, which is defined 
by 


c k (t G^ “ ir a l W. 


K i K j_ 


T . T . 
si sj 


W G i (y ik ) G i ( yjk ) 


( 11 . 2 ) 


vhere and Kj are the sensitivities of the two antennas in 
ielvins/jansky . This formulation of the amplitude is based on 
the assumption that the correlated flux and antenna sensitivity 
are essentially constant across the (2-MHz) bandpass. The factor 
bjj accounts for the two possible normalization options and is 
given by 


bN = N-j- for no normalization (11.3) 


= 1 for N t normalization (11.4) 


In Eq. 11.2, bjj has been made a function of tg to account for 
the fact N-£ might change from interval to interval as a result 
of lost data. 

For the special case of a flat power spectrum Sp and identical 
bandpass shapes at the two stations, the "bandpass factor" in 
Eq. 11.2 can be converted through the use of Eqs. 6.12 and 8.12 
to the form 
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(11.5) 


W G i ( rik> G i^ik ) s 6 <V 


D, 


G 2 


where 


G 2 



G 2 (o))dco 


( 11 . 6 ) 


in which doppler shifting has been neglected. For a rectangular 
bandpass one can readily show that the bandpass ratio in Eq. 

11.5 is 1.0 so that, in this approximation, the fringes are 
given by 


< h c ( "k' t G ) > 


S - 
C TT 


a l b N 



i »c ( "k' t G ) 


(11.7) 


This crude form will be useful in Section XIV in defining 
default values for amplitude in the fringe model. 
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SECTION XII 


CONVERSION TO COMMON GREENWICH TIME 


To make all components of delay a function of "common Greenwich 
time," tg, this section consolidates the time tags that 
characterize fringe phase. The (true) time tags tpi and tpj 
in Eq. 10.6 are related to common Greenwich time tg through 
Eqs. 5.6 and 6.4, which will yield the relation 


t pk - t G “ T bk “ T Ik (12.1) 

To convert the time tags of all quantities to tg, expand x^i and 
X(jj in Eq. 10.6 about tg to obtain 

AT R (t G ) = Axg(tg) + Ax t (tg) + Ax c (tg) - x u + T a + Ax cs (12.2) 

in which Eqs. 6.1 and 6.20 have been used to establish the 
definitions for the residual geometric delay: 

V [Vv - 5 i (t S > ] V [^nj (t O> - ^ni (t G > ] 

*yv 1 c 2 5 (12 ' 3> 

for the residual troposphere delay: 


AT t <t G ) = 


x tj (t G ) " 


T ti (t G ) “ 


m / 1 \ 
T tj ^ t G 


m x 


(12.4) 


and for the residual clock offset: 


4T c <t G ) = x cj (t G ) " t ci (t G ) ' T cj (t G> 


, m /. v 

+ t ci (t G ) 


(12.5) 


In Eq. 12.5, the 3-lag shift in Aj discussed in Section VIII 
has been taken into account so that Aj - A^ can be set 
equal to zero. The term Ax cs arises in the expansion about 
tg and is given by 
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( 12 . 6 ) 


Al. 3 X 
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«jj<V 
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x ij > ■ 


T di ( V 


(x 


bi 


T Ii } 


Thus, to apply a correction to remove this last term, one must 
determine the instrumental delays xjk and the absolute bit-time 
errors x^k that give the offset of the sampler clocks from true 
time. Even if these quantities sum to 10 nsec, the correction 
would be less than 2 psec/sec • 10 nsec = 20 psec . After such a 
correction, the only clock synchronization effect remaining in 
the observed delay will be the additive term Ax c , the offset 
between phase-cal clocks at the two stations. If these time-tag 
corrections are carried out, the time tag of the observables can 
be precisely treated in subsequent modeling as though it referred 
to a true time of t q at both stations. The correction would 
probably be applied to the overall observables for each observa- 
tion (1-10 minutes), and would be computed on the basis of aver- 
age delay rates for the observation (see Section XIV for a 
discussion of the extraction of final observables). 

The fringes have been corrected for instrumental and processing 
effects and are theoretically given by Eq. 11. 1, in which fringe 
phase y c is given by Eq. 10.5 with residual delay Axr given by 
Eq. 12.2. They will now be transformed from the RF-f requency/ 
time domain to the BWS-delay/f r inge-f requency domain where signal 
detection and observable extraction can be readily carried out. 
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SECTION XIII 


TRANSFORMING TO THE BWS-DELAY AND FRINGE-FREQUENCY DOMAINS 


When an accurate correlator model is used, the counter-rotated 
fringe phase in Eq. 10.5 is nearly a linear function of both time 
and frequency. Thus, if the frequency slope (i.e., residual BWS 
delay) and time slope (i.e., residual fringe frequency) were 
relatively accurately known, the fringes could be counter- 
rotated over both dimensions with a linear phase model and 
coherently summed over time and frequency. Those slopes, 
however, are generally not known accurately enough to perform a 
coherent sum. Thus, counter-rotation for a range of "trial" 
slopes is carried out by means of an FFT over frequency and an 
FFT over time. These FFTs transform the fringes from the 
RF-f requency/time domain to the BWS-delay/ fringe-frequency 
domain. For output points (slopes) near the actual BWS delay 
and fringe-frequency, the fringes will coherently sum to form a 
peak in output amplitude that will allow detection and observable 
extraction. Appendix D and Figures 5a and 5b outline the trans- 
formations in terms of data structure. 


To carry out these transformations, the fringes in Eq. 11.1 are 
counter-rotated by two consecutive FFTs as follows: 


R’ (x f ,co f ) 
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(13.1) 


where, for standard fast Fourier transforms, the time and 
frequency variables are both zero for the first point: 



(13.2) 


t G t G t l (13.3) 

where co^ and t^ are the values for the first frequency and time, 
respectively. 


TRie sum over frequency is normally performed over the bins from 
all channels in a band (e.g. , S or X) , while the sum interval 
over time is selected by the user. For a standard FFT, the 
"trial" values for (BWS) delay Xf and fringe frequency cof 
are given by 
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where Axf (Acof) is the separation between the delay (fringe- 
frequency) points. (Negative "trial" values of BWS-delay and 
fringe frequency reside in the upper half of a standard FFT 


13-1 



output, as noted in Appendix D.) As illustrated in Appendix D, 
the separation for the BWS-delay coordinate is of approximately 
1/4 to 1/2 the inverse of the RF spanned bandwidth while the bin 
width in the frequency direction is approximately 1/4 to 1/2 
the inverse of the integration time. 

Because the bin spacing is 0.5 MHz in RF frequency, the FFT 
output is cyclic as a function of BWS delay, with a cycle 
interval of 2 nsec. Thus, to obtain absolute delay measurements, 
the a priori error on the calibrated delay must be no greater 
than 1 psec. Typically, a priori errors due to geometric, 
tropospheric and ionospheric delays will not be large enough to 
prevent resolution of this n«2 nsec ambiguity. The offset 
between station clocks, however, can potentially be a large 
unknown value. In this situation, a relatively accurate esti- 
mate of residual delay can be obtained from the single-channel 
delay obtained directly from the lag-domain fringes generated by 
a strong source. This offset, however, must first be corrected 
for the instrumental offsets t j and x^. With stable clocks, 
analysis of a single observation in this manner should usually be 
adequate to determine the absolute clock offset and thereby allow 
removal of the n»2 nsec ambiguities for a whole session. 

As illustrated in Figs. 5a and 5b, the FFT over RF frequency is 
carried out first. To reduce memory and computation require- 
ments, only a subset of resulting BWS delay points, selected on 
the basis of the a priori delay, are retained to be transformed 
over time (see Appendix D) . 

To compensate for the two origin offsets in Eq. 13.2 and 13.3, 
the output of the FFTs is counter-rotated at each point as 
follows: 


R(x f ,to f ) 


i[ (co z -a) 1 )T f +(t z -t 1 )a) f ] 


R' (x f ,io f ) 


(13.6) 


This counter-rotation places the "origin" of the FFT frequency 
variable at a reference frequency co z and the "origin" of the FFT 
time variable at the reference time t z . Both reference values 
can be specified by the user but should be computed according to 
the expressions in Appendix C. Both values are analytically 
enforced by the counter-rotation in Eq. 13.6 and by the least- 
squares model explained in the next section. Based on Eqs. 

10.5, 11.1, 13.1, and 13.6, the expectation value of the 
resulting fringes is given by 


<R(x f ,u) f )>= 


k,t. 


s c c k e 


H W W < ‘V’z > x £- ' ! V* z > ' “£. 


(13.7) 


where Ax R is given by Eq. 12.2, and c^ is given by Eq. 11.2. 
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This fringe function peaks in amplitude at the values for xf and 
wf that are closest to the actual BWS delay and delay rate, as 
illustrated in the next section. For stable instrumentation, the 
shape of the fringe function is a sinx/x form along the fringe- 
frequency axis. Along the BWS-delay axis, the fringe function 
has a main lobe centered on the true delay with a sidelobe struc- 
ture dependent on the placement of channels within the spanned 
bandwidth, relative amplitude of the channel fringes, and the 
coherence between channels. 


For the transformation approach outlined above, amplitude can be 
lost if a priori values for baseline and source position are 
very inaccurate. Since the residual fringe frequency in the 
passed fringes scales with channel frequency, each channel will 
have a different fringe frequency. Thus, in the double FFT over 
time and RF frequency, the output of the time FFT will not be 
perfectly coherent when summed across channels. The smaller the 
residual delay rate, the smaller the amplitude loss due to this 
effect will be. A limit can be placed on such amplitude loss by 
adjusting the integration time (Tj) so that the change of 
residual fringe frequency across channels is considerably less 
than the width (2/Tj) of the main lobe of the FFT's frequency 
filter. As an example, the following approximate formula for 
the maximum integration time applies when the BWS channels are 
uniformly spaced across the RF band: 


2.5 


V^A 


max 


(13.8) 


4f s iT * 


where is the maximum fractional decrease (e.g., 0.01) in 

amplitude to be allowed, Af g is the spanned bandwidth, and Ax^ 

is the maximum expected residual delay rate. Other distributions 
of channels across the RF band may have significantly different 
values for maximum integration time. 


To illustrate the limits set by Eq. 13.8, let the desired inte- 
gration time be 1000 sec, the approximate limit set by the 
instability of the current frequency standards. For a spanned 
bandwidth of 400 MHz and a maximum amplitude loss of 10%, Eq. 

13.8 would limit the residual delay rate to about 2 psec/sec or 
less. A delay rate error of 2 psec/sec corresponds to a baseline 
error of about 8 m and a source position error of about 0.2 arc- 
seconds (for a 9000-km baseline). These loose requirements on 
the a priori delay will almost always be satisfied if the 
baseline and source have been previously measured by the VLBI 
technique. In the occasional cases for which the source or 
baseline is being accurately measured for the first time, 
preliminary "calibration" observations may be necessary (e.g., 
by "packing" the channels side-by-side for a few observations or 
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observing the source first on a shorter baseline) . As VLBI 
matures, a smaller and smaller fraction of sources and baselines 
will have poor a priori delays. 

These accuracy requirements for a priori baseline and source 
positions are based on a uniform distribution of channels across 
the spanned bandwidth. As another example, one can use the 
channel frequency values listed in Appendix E, which are closer 
to the distributions used to optimize delay measurements. For 
this distribution, the limits on a priori that are presented 
above must be tightened by about a factor of 1.3. 

The amplitude loss in Tensor discussed above occurs when the a 
priori delay model is very poor. Since special methods can be 
employed to treat such cases, the only impact of this limitation 
is the extra effort that might be required the first time an 
inaccurate baseline or source position is used. Plans have been 
formulated for correcting this limitation in Tensor, thereby 
increasing the allowed range of uncertainty for the a priori 
delay model. 
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SECTION XIV 


OBSERVABLE EXTRACTION 


The first step in extracting information from the fringes 
produced by the operation in Eg. 13.6 is to find the maximum 
amplitude (RSS of real and imaginary parts) as a function of BWS 
delay (xf) and fringe freguency (tof ) . A priori information can 
be used to restrict the range over which the search for the 
maximum is carried out. Since the fringes that are input to the 
FFT have a frequency spacing of 0.5 MHz, the output fringes will 
repeat every 2 jasec as a function of BWS delay. The a priori 
information should be accurate enough to prevent this n*2 nsec 
ambiguity. Even if the a priori error in delay is too large 
(> 1 nsec) due to clock error, the uncertainty in the other 
components in the delay will usually be accurate enough to 
select a common (but arbitrary) ambiguity for all observations. 

In subsequent multiparameter fits, this compromise will not 
degrade estimates of geophysical and astrometric parameters, 
but will, of course, prevent synchronization of station clocks. 

The FFT output point at which the maximum amplitude occurs will 
be called the peak point. Data in the vicinity of the peak point 
are fit by least-squares to extract the fringe observables. Only 
fringe values within 1 to 3 intervals of the peak point (9 to 49 
points in all) are included in the fit. The model for the 
fringes in the least-squares analysis is given by 


i<t> 


_ , * , ^ s S i (w^-to,) (x o -x^) + (t_-t _)(w.T 

V T f'“f> =A s e f- 1 

k,z g 


( 




(14.1) 


vhere t z , co z are the reference time and frequency used in 
Eq. 13.6. The four solve-for parameters are # amplitude A s , BWS 
delay x s , RF phase <j> s , and phase-delay rate x<t> s . 


Hodel values for the relative amplitudes c m fc can currently 
3>e specified in two ways. In the default mode, they are 
automatically set equal to a crude value that is valid for a 
system in which all channels are identical (e.g., have the same 
correlated flux, the same rectangular 2-MHz-wide bandpass 
function, the same antenna sensitivity, and the same system 
temperature) . That default value is given by 



(14.2) 
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As one can see by comparing the default value with the amplitude 
in Eq. 11.7, the default amplitudes have been selected so that 
the estimated fringe amplitude in the default mode will be 
approximately given by 
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(14.3) 


where the temperatures are the common temperature for all 
channels. When- accurate delay and amplitude measurements are 
the goal, one has the second, more accurate option of supplying 
to the post-Tensor software numerical values for the relative 
amplitudes. With this option the user has the maximum flexibil- 
ity in modeling the fringe amplitudes, regardless of his approach 
for obtaining them. One approach is to measure the quantities 
necessary to evaluate the theoretical expression for amplitude in 
Eq. 11.1. For example, if the correlated flux is the same for 
all channels across the spanned bandwidth, the relative ampli- 
tudes are theoretically given by Eq. 11.2, which requires esti- 
mates for bandpass shapes, system temperatures and antenna sensi- 
tivities for each channel. (In very accurate applications, it 
may be necessary to upgrade Eq. 11.2 to account for intrachannel 
errors in relative bin amplitude that were introduced by the 
limited-lag transform but neglected after Section VIII.) Under 
these assumptions, the estimated amplitude will be equal to the 
correlated flux (As = S c ) . An alternate, less theoretical 
approach for estimating relative bin amplitudes is to separately 
extract the bin amplitudes for a strong source with constant 
correlated flux across the spanned bandwidth, and to use their 
relative magnitudes along with a scaling factor to compute the 
values for c m jj. An error in the scaling factor will cause a 
scale error in the estimated fringe amplitude A s but will not 
degrade the estimate of BWS delay, delay rate, and phase. If 
the correlated flux of the strong source is known, the scaling 
factor can also be obtained from the calibration run. 

If accurate values for relative bin amplitude are supplied and 
if the observed fringes are sufficiently coherent over frequency 
channels, one can show (Appendix C) that a best fit will be 
obtained when the model phase (phase in Eq. 14.1) is adjusted so 
that it approximates (in a least-square sense) the actual phase 
(the phase in Eq. 13.7): 


<“k - “z ><T s “ x £ ) + (t G - V (l V<t>s “ "£> + +s = 

“k iT R <t G ) + VVV + V'VV - <“k - V T f 
- (tg - t z ) co £ + 2nir (14.4) 
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where 2nir represents the integer-cycle ambiguity in fringe phase 
and where AxR(tg) is given by Eq. 12.2. If one cancels like 
terms on either side and expands the structure phase 4 >b about 
(w z , t z ) and Atr about t z , this equation can be rewritten as 


♦s + <to k - “z> x s + (t G - V “k x 4>s “ 

VW + <“k - “z> ix R ( V + (t G " V “k Ax R (t z> 

8<J> B (t0 z ,t z ) 8< * > B^ ( °Z ,t Z > 

+ V'W + ( “k - "z> HT + (t G ~ V It 


k e W 
to,. 


(14.5) 


where the last term is the ionosphere phase 4> e expressed in 
terms of frequency and differential electron content AN e (t). 
Theoretical expressions for the three estimated parameters can 
be obtained by substituting the right-hand-side of this equation 
(which is the expansion of Eq. C.14 in Appendix C) in place of 
<j; c in Eqs. C.19, C.20 and C.21, which yields the following 
approximate expressions for BWS delay, phase-delay rate and 
phase-delay: 
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where co T , co^, and co<j> are effective frequencies to be used for 

for computation of the ionosphere effect and AN e and AN e are 
computed on the basis of electron content. Expressions for AN e , 

AN e and the three effective frequencies are derived in Appendix E. 
Residual delay Atr is given theoretically by Eq. 12.2, but evalu- 
ated at reference time t z . 

In most applications, the model fringes in Eq. 14.1 should be 
carefully constructed. When accurate estimates of BWS delay are 
required, for example, the relative magnitudes of the bin am- 
plitudes c m k should be accurately specified so that the shape 
of the "delay resolution function" closely matches the actual 
shape. The use of the accurate relative bin amplitude also 
simplifies the problem of computing effective frequencies for 
dual-band ionosphere calibrations, as explained in Appendix E. 

Another important consideration in assessing the accuracy of the 
fringe model is the phase coherence of the observed fringes. 

Since phase coherence was an important goal when the station 
instrumentation was designed and developed, the fringe model in 
Eq. 14.1 has been based on the assumption that the instrumental 
contribution to fringe phase is coherent (i.e., perfectly linear) 
as a function of time and RF frequency. It is likely, however, 
that the actual level of coherence will not be adequate in some 
applications. Given a strong source, a test of coherence can be 
carried out by extracting the fringe phase at each time point 
(e.g., each second) and for each frequency bin for each channel 
and then fitting all of those phase values with a phase model 
that is linear in time and frequency (e.g., as in Eq. 14.1). 

Plots of the residuals will reveal nonlinear trends as a function 
of time and frequency. Incoherence over time usually cannot be 
modeled and corrected since it is typically a consequence of 
instability in the frequency standards. Incoherence over 
frequency, on the other hand, is likely to consist of the same 
nonlinear pattern from source to source. If that is the case, 
bin-specific phase corrections can be computed on the basis of 
the aforementioned "coherence residuals" across frequency and 
added to the phase of the model fringes. Such a correction 
would make the shape of the model "delay resolution function" 
more accurately represent the actual shape and would lead to 
more accurate estimates of BWS delay. 

An example of a fit to Tensor output is shown in Figures 7a and 
7b. These data, which were reduced by T.G. Lockhart using his 
post-Tensor software, were obtained on an OVRO/Mojave baseline 
for the source 3C84 . Before the Tensor FFT's, the data in the 
time and channel-frequency domains consisted of 30 time points 
spaced by 2 seconds and 8 upper-sideband channels with L.O. fre- 
quencies of +( 0 , 10, 40, 100, 210, 290, 340, 360) MHz relative 
to 8210.99 MHz. A 2048-point FFT over channel frequency, 
followed by a 64-point FFT over time, produced the data points 
shown in Figures 7a and 7b. To reduce data volume, only two 
slices of the full two dimensions are presented. The slices 
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pass through peak amplitude along the delay axis and along the 
delay rate axis. Note that the amplitude has been normalized to 
correspond with the "units" of estimated amplitude, A s , in 
Eq. 14.3, with b n = 1 . The solid curves in the figures repre- 
sent the fit function given in Eq. 14.1. A reference frequency 
of f z = 8380.74 MHz was used in the fit. The sharp 1/2 cycle 
changes in phase vs. delay and delay rate are due to abrupt 
changes in quadrants when the amplitude becomes very small. The 
fit as a function of delay rate is very good for both amplitude 
and phase. On the other hand, the fit as a function of delay is 
fairly good in the main lobe, but deviates considerably in the 
sidelobes. These deviations, which should be eliminated by 
improvements in hardware and/or better modeling, are probably 
due to non-linear variations of fringe phase across frequency 
channels or by unmodeled variations in fringe amplitude across 
frequency channels. 
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SECTION XV 


MODEL RESTORATION 


After the (residual) observables have been extracted by the 
least-squares fit, the post-Tensor software computes a model 
delay (or delay rate) based on the processor model and adds this 
model delay (or delay rate) to the estimated observables in 
order to form total observables free of processor modeling. 

Two approaches were considered for computing the restored model: 
analytical recomputation of the correlator model for each delay 
component (i.e., geometric, tropospheric, and clock) or inter- 
polation between the actual numerical delay values used by the 
correlator. Since analytical recomputation of the processor 
models is prone to error and inconsistency, the numerical 
approach was selected. The correlator passes for each station 
an array of values for the total delays x m i, x m j used to calcu- 
late the quantized bitstream offsets. These array values, which 
the correlator computed on the basis of Eq. 6.2, are separated 
by 2-4 seconds. The post-Tensor software interpolates between 
the passed values to obtain model values at desired times. 

Were it not for the "retarded baseline" ("transit time") compli- 
cation in the geometric component, the interpolated model delay 
values could simply be differenced between stations and added to 
the residual delay. For geophysical/astrometric applications, 
however, the post-Tensor software must restore the model in a 
form that will lead to a geometric component that is based on the 
definition commonly assumed in subsequent multi-parameter fitting 
software: a wavefront arriving at antenna i at true time t will 
reach antenna j at true time t + Xg(t), where Xg(t) is defined to 
be the geometric delay at true time t. This definition of the 
geometric delay incorporates the effect of the motion of station 
3 during wave transit, an effect that is not present in the 
straight difference between stations. 

The preceding definition can easily be converted to an equation 
for computing geometric delay in terms of station locations and 
source direction. The ph^se of the incoming signal is given in 
terms of the w^ve vector^k and frequency oj by the exponent in 
Eq. 5.1. Let x^(t) and xj(t) represent the positions of the two 
antennas as a function of J true time. The phase arriving at 
antenna i at time t will reach antenna j at time t + Xg when 
the following phase relation is satisfied: 


Ic • x.(t) - cot = lc • x.(t + x ) - co(t + x ) (15.1) 

-*■ j y y 
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Rearrangement of terms produces an expression for Tg: 

x g (t) = | • [xj(t + x g ) - x i (t)J (15.2) 


where k is the unit vector in the direction of k. The difference 
in station locations in brackets is equal to the baseline vector 
between stations except for the small offset by Xg in the argu- 
ment of xj . This argument offset is the origin of the term 
"retarded baseline" correction. 

To illustrate the magnitude of the retarded-baseline correction, 
expand x j ( t ) about t : 


x.(t + x g ) * Xj (t) + Xj x g 


(15.3) 


in which higher order terms have been neglected. For geometric 
delays less than 0.03 sec, the higher order terms contribute 
less than 



i | 2 

|x j |x g 


(15.4) 


< 0.1 psec 


(15.5) 


which is negligibly small. Thus, the magnitude of the retarded- 
baseline correction is bounded by 


iT g S c |x g (15 ' 6) 


< 3 X 10 6 * 0.03 sec (15.7) 


< 90 nsec 


(15.8) 



for intercontinental baselines. Since this effect is relatively 
small, relatively inaccurate estimates for source positions and 
baseline vector provide accurate estimates for this correction. 
For example, a 100-m error in baseline vector or a 2 arcsecond 
error in source position will result in an error of less than 
0.1 cm in the correction. Thus, the retarded-baseline effect 
can be accurately introduced in the final observable so that the 
conventional definition of geometric delay will apply to the 
geometric component of the total delay. 

The expression for geometric delay can be cast into a more prac- 
tical form by replacing Tg on the right-hand side of Eq. 15.3 
by the symmetric geometric "delay" 


Tg(t) = | • [xj(t) “ ^(tjJ (15.9) 


Since the retarded baseline effect is less than 30 m, this sub- 
stitution results in errors less than 0.2 psec . When Eq. 15.3 
is combined with Eq. 15.2, the geometric delay becomes 


T g (t) ~ T g (t) + Ax g (t) 


(15.10) 


where the symmetric geometric "delay" x£j is given by Eq. 15.9 
and the transit-time correction by 


T*(t) 



(15.11) 


In order to produce observables satisfactory for both source 
structure and geophysical/astrometric applications, the model 
delay is restored in two steps. The first step restores a model 
delay whose geometric component is of the form of Eq. 15.9. This 
leads to total observables that are symmetric with respect to 
interchange of stations, as required by source structure work. 
The second step adds a transit-time correction of the form of 
Eq. 15.11, which produces total observables based on the conven- 
tional definition of geometric delay, as assumed in geophysical/ 
astrometric applications. 

In the first step, the restored model is computed as 


WV ■ *mi<V - *mj<V 
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in which the unretarded model delays are obtained by interpola- 
ting between the numerical values passed by the correlator. 
Total observables, essentially free of correlator modeling, are 
obtained by adding the numerically interpolated models to the 
residual observables produced by least-squares analysis: 
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where x^ is the numerically computed time derivative of x^. 

One can obtain a theoretical expression for the total observables 
of step 1 by combining the theoretical forms for the observables 
in Eq. 15.13 to 15.15 with the theoretical expression for the 
residual delay obtained from Eqs. 12.2 to 12.5 and Eqs . 14.6 
to 14.9: 
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in which xgj is the actual (symmetric) geometric delay given by 

Eq. 15.9 with k = k a , and where and x c are the differences 
between stations of the total actual troposphere and clock terms, 
respectively, in the order j minus i. These total observables 
are passed to source-structure software for further processing. 

The second step in model restoration produces observables suit- 
able for geophysical/astrometric applications. In this step, a 
model delay that accurately approximates the transit-time correc- 
tion in Eq. 15.11 is added to the observables computed in step 1 
(Eqs. 15.13, 15.14, and 15.15). As indicated earlier in this 
section, a relatively inaccurate model for source position and 
station location will not significantly degrade final delay accu- 
racy. The resulting observables, which will be denoted by a 
superscript^ 2 ) , are still given theoretically by Eqs. 15.16, 
15.17, and 15.18, except that x R (D becomes 
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where xg(t z ) is now the actual geometric delay based on the 
conventional definition, as in Eqs. 15.2 or 15.10. 

As shown by Eqs. 15.16, 15.17 and 15.18, the total observables 
are a sum of terms due to geometry, troposphere, clock, charged 
particles, source structure, and three "instrumental" effects. 
When absolute clock synchronization is a goal, postprocessor 
software will apply corrections to remove the phase-cal cable 
delay x u and the "antenna" delay xa- In future ultraprecise 
applications, the time-tag term Ax cs and the source-structure 
terms involving <$>g might be removed also. After these correc- 
tions have all been applied, the final three observables (BWS 
delay, phase-delay rate, and phase delay) are each a sum of terms 
due to geometry, troposphere, charged particles, and clock off- 
sets, where all terms can be modeled as though they were observed 
at true time t z . As mentioned in Section V, the geometric 
delay is modeled relative to the effective positions (inter- 
section of axes) of the two antennas. This set of observables 
is extracted for each RF band. 
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SECTION XVI 


SUMMARY 


This report has presented radio interferometry theory as it 
applies to the BLKII processor, including many detailed instru- 
mental and processing effects that are important in source-struc- 
ture and geophysical/astrometric measurements. The analysis has 
covered all major data reduction steps from signal reception to 
the output observables of amplitude, BWS delay, phase delay and 
phase-delay rate. Even though the BLKII processor is vastly 
superior to its predecessor in a number of ways, a number of 
additional improvements are being considered, as noted in the 
introduction. We plan to implement many of these improvements 
as resources permit. 
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APPENDIX A 


A COMPLETE DERIVATION OF THE CROSS-CORRELATION FUNCTION 


This appendix presents a complete derivation of the cross- 
correlation function, taking into account the statistics and the 
spatial distribution of the emissions of the natural source. 

The radio noise generated by a very distant extended natural 
source can be expressed as a superposition of plane waves in the 
form 


E(x, t ) 


n 00 - -> 

A( k , w) e i “ [t -' t,x/o! dudfl + c.c. 

K 0 


( A . 1 ) 


where E(x,t) is the electric field at position x and true time t. 
For simplicity, the electrig £ield is assumed to be linearly 
polarized. The function A(k,x) is the Fourier ^amplitude at fre- 
quency co for the wave received from direction k, while c.c. 
denotes the complex conjugate. As in Section V, all quantities 
are measured with respect to a geocentric frame. The wave 

direction k can be expressed as a function of two parameters 
and, in the case of right ascension (r.a.) and declination, 
becomes 


k = - (cosS cosa, cosS sina, sinS) (A. 2) 

where a, S are the apparent right ascension and declination 
relative to the true equatorial coordinates of date. The quan- 
tity dfl represents a differential solid angle such as cosS da dS 
in the case of right ascension and declination. In the following 
steps, we will use the wave vector k = oak/c. 

The electric field detected at antenna j is given by 

r r i[ajt-£«x. (t) ] 

E . ( t ) = E(x . ( t ) , t ) = J / A(k,co) e J dajdn + c.c. 

3 3 j i 0 (A. 3) 

where xj(t) is the position of antenna j as a function of true 
time. {.See Section V for a definition of antenna location.) 

This signal will be received by the antenna and passed through 
various filters and mixers. The modeling of the composite 
effect of instrumentation is discussed in Section V. Following 
that derivation, the signal component due to the natural source 
can be written in the form 
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+ c . c . ( A . 4 ) 

where the doppler-shifted passband argument is given by 

(A. 5) 


yj = to (1 - k • xj/c) 


and where instrumental phase terms are given by 


a Ij<“' t pj > = “hj t pj + *hj + “ T aj + 4, Ij <y j > <A ' 6) 


As in the text, tpj is the true time that corresponds to bit 
time tbj ■ (See Eq. 5.6.) 

We will first compute the expectation of the product of two 
analog signals as described by Eq. A. 4 and then compute the 
normalization. The expectation value of the voltage product 
is given by 


<vf(t M ) v?(t bj » - 


I FI I 

■ J n J \r J n 


oo 


A(k,to)A*(k 1 ,co 1 ) G^y^ Gj(yj) e 1 ’*' dco dfl dco 1 dQ’ 


i ft r 


+ 1111 A(k,co) A(k',co') G^(y^) G j(Yj) dco dfl dco' dfl' 

k‘ 0 k 0 


+ c . c . 
where 


(A. 7) 


yj = co' (1 - k' • Xj/c) 


(A. 8) 


v = »t pl - t • 5i<i pi ) - “'tpj + *' • Xj <t pj ) 


- a n <-,t p .) + <*!}<<■>'. t pj > 


(A. 9) 
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and Yb is a similar expression that will not be needed. 

The quantities t p k are the offset values for the true time 
corresponding to t^k' as in Eqs. 6.8 and 6.9 in the text 
(i.e., tfcj = tfcj + A implies t p j = t p j + A). 

We will assume that the natural source is completely incoherent, 
which means 


<A(k,io) A*(k' , to ' )> 


1 

2 S p (k,co) 


5 (k - k ' ) 6 (to - u 1 ) 


(A. 10) 


where S p (k, to) is the power spectrum for direction k and 5(z) 
is a Dirac delta function.* That is, emissions from different 
areas of the source are uncorrelated. Furthermore, the noise 
emitted by a given area is stationary and therefore possesses 
uncorrelated frequency components. Since the signal is real, 
we know 


A(k,co) = A*(k, - co) 


(A. 11) 


This relation and Eq. A. 10 give 


<A(k,co) A(k',w , )> = \ S (k,co) 5 (k - k 1 ) 6 ( oj + co') (A. 12) 


= 0 for co and co' > 0 


(A. 13) 


Under these assumptions, Eq. A. 7 becomes 


r r 00 

<v f(t bi ) V^(t bj » = \ I J 8 <i.«) G^y.) e I <toKin 


k o 
+ c.c 


(A. 14) 


where 


y n = to (1 - k • x n /c) n = i or j 


(A. 15) 


*For two particular direction parameters (13, y), the 6(k - k') 
delta function denotes 6(13 - 13') 6(y - y'). Furthermore, we 
will require 13 and y to satisfy the relation d!3 dy = dD. 
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Yj = w (t pi - t p j) + t • S r - « ri («.t i> + (A. 16) 


The A(k,aj) A(k',oo) term has dropped out since w = — co ‘ is not 
covered in the region of the integration. We have also defined 
a "retarded baseline": 


K = s j (i P j ) - 


(A. 17) 


Let two particular parameters (8, y) define the direction vector. 
Suppose the brightness distribution is very narrow about some 
reference direction k a : 


k = k <8, 


V 


(A. 18) 


If the brightness distribution is sufficiently compact, we can 
approximate the wave direction by 


k = k + — 
K K a + 813 


<I3-B a ) + |f 


(Y-Y a ) 


(A. 19) 


where the partials* are evaluated at the reference point 8 a , y a . 
We then obtain 


r™ _ _ i<j) 

< v ? ( t bi ) V?(t b .)> = \ J R(u,v,to) G i (y i ) G.(y.) e f dco + c. 


c . 
(A. 20) 


where 


y n = “ (1 ‘ k a • V c) 


n = i or j 


(A. 21) 


*f = 03 ( *pi ~ t pj ) + * a * *r “ a Ii (w 'Si ) + a Ij (a) '^pj ) (A - 22) 


^Structure coordinates are usually defined as B = (a - a c ) cos 6 C 
and y = S - S c where (a c , 6 C ) is the origin of structure coordi- 
nates relative to r.a. and declination. With this definition, 
the 8 partial and y partial are in the direction of increas- 
ing right ascension and declination, respectively. 
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and where 


R B ( u , v , to ) 



Sp(f3, y,w) 


-2iri[u(I3-fi _) +v(y-y )] 
e a a dfldY 

(A. 23) 


for which 


u = 




(A. 24) 





(A. 25) 


X = 2irc/oj 


(A. 26) 


With this definition, u and v are the components of the baseline 
vector projected on the plane of the sky. 

Two approximations have been made in Eq. A. 23. First, the weak 
k dependence in the bandpass functions has been neglected so that 

the yi, yj values have been evaluated at k a . Second, the 
limits of J the (13, y) integration have been extended to infinity 
under the assumption that the brightness distribution is very 
compact and terminates the integration. The function R(u,v,to) 
is referred to as the brightness transform. A more detailed 
treatment of the brightness transform and of structure effects 
is given in Ref. 7. 

The fringe visibility is the normalized modulus of the brightness 
transform and is a measure of the spatial extent of the source 
relative to the resolution of the interferometer for the given 
observation. It is equal to the ratio of "correlated flux" to 
total flux and is defined by 


v f (u, v,co) 


Rg ( u , v , to ) 
R B ( 0 , 0 , to ) 


(A. 27) 


R 0 (u, V, to) 

y u >"' 


(A. 28) 
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where Sp(co) is the total power emitted by all parts of the source 
at frequency to. For a point source, the fringe visibility is 
unity for all baselines. For a diffuse source, the visibility 
is unity for a "null" baseline and decreases on the average as 
the baseline length increases. The decrease can become 
substantial for baselines satisfying 


-*■ 

B 


> 


X 


~ 2d. 


(A. 29) 


where d s is the source diameter in radians and X is the radio 
wavelength . 

In addition to this normalization of the brightness transform, 
we must normalize the voltage signal in Eq. A. 4. Using Eqs . A. 10 
and A. 13 one can easily show that the expectation of the squared 
voltage is given by 


_ 00 

« V n , 2 >=i-I V*' u) G n (y n> 


dudU 


(A. 30) 


k ^o 


-f 


Sp(u) G^(y n ) do) n = i or j 


(A. 31) 


where we have again neglected the weak dependence of y n on k and 
evaluated y n at the reference position k a . Using Eqs. A. 6, A. 17, 
A. 20, A. 22, A. 28, and A. 31, we can now compute the cross- 
correlation function for the normalized voltage components for 
the natural source: 


<9i<i b l) = 


< V i (t bi - \u - tb > v ?< t bi - W> 

T rr 1 rt 2 -1 I / 2 

[< (V?) > <(V?) >] 


(A. 32) 


Vf r = i4> f 

/ S p ( M ) G.(y.) G j (y j ) e do) + c.c. 

IN J ry 


(A. 33) 
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where 
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= f S (cu) G?(y. ) doj f 

*/q P 11 


00 


S p (co) Gj (yj ) do> 


(A. 34) 


and where 


*f m ♦h “ “hj \ij + + lb) + *i + 


+ k a * Cx j {t pj “ W - x i (t pi - T mi - 

+ »(T b + Tj + x a - x mi + - lb) 


(A. 35) 


in which x^, xj, x a , and <J>j are differences between station 
quantities in the order j minus i. 

As explained in the text in Eq. 6.17, the timetag difference, 
tpi ~ tpj , is a combination of the difference in the errors in 
the bit times of the two stations and the difference in instru- 
mental delays. We have pulled vf out of the integral, recog- 
nizing that resolution of the source will vary very little across 
a 2-MHz passband located at RF. A phase term 4>g has been added 
to the fringe phase 4>f to represent the phase term from the 
brightness transform R. Further, we have used the definition of 
the instrumental phase from Eq. A. 6, properly offset in time as 
indicated in Eqs . 6.8 and 6.9. Note that this expression for the 
fringe phase is the same as the expression (Eq. 6.13) assumed 
without proof in the text, except for the phase shift <t>e due to 
charged particles. The charged-particle phase shift has been 
added in an ad hoc fashion in the text to allow for such effects. 

As indicated in Eqs. 6.5 and 6.7, the final expression (Eq. 6.10) 
for the cross-correlation function of an extended source is 
obtained by multiplying Eq. A. 33 by 2/ir and by a noise tem- 
perature factor. One can easily show that this result for an 
extended source is the same as that for a point source except 
that (a) an amplitude factor (vf < 1) must be included, (b) a 
phase shift due to structure must be added to the phase, and 
(c) the wave direction is set equal to the reference direction 
k a . 
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APPENDIX B 


THEORY FOR CALIBRATION-TONE PHASE 


This appendix treats the subject of phase calibration with a 
tone generator and derives an expression for the stopped tone 
phase. 

The calibrator signal, which consists of a sum of constant- 
freguency tones uniformly spread across the passband, is injected 
at a low power level near the front of the instrumentation and 
is imbedded in the radio source data. At the correlator, each 
tone is counter-rotated to nearly zero frequency by digitally 
mixing each single-station bitstream with a quadrature sinusoid 
whose frequency closely approximates the baseband frequency for 
that tone. After the resulting tone phase has been extracted 
and differenced between stations, it is subtracted from the 
interferometer (fringe) phase. Since the tones are subjected 
to the same instrumental effects as the natural source signal, 
instrumental effects between the injection point and the recorder 
cancel in the final difference. (This report will not consider 
difficulties encountered in calibration of instrumental effects 
before the injection point.) 

Two types of phase calibrator are currently in use — the original 
Rogers version and the later Sigman version. The Rogers phase 
calibrator generates tones with a frequency spacing of 1 MHz. 

Its cable stabilizer records, but does not simultaneously 
remove, variations in the electrical length of the "uplink" 
cable. An average of the recorded variations is applied in 
subsequent processing to correct the composite observables 
obtained for an observation (e.g., over 3 minutes). With the 
Sigman calibrator, the frequency spacing of the tones can be set 
by the user to one of the values given by 5/N MHz, where N = 1, 

2, ... or 99. The 5/N-MHz selection of zero crossings is 

automatically related to the station 1 pps signal in order to 
facilitate clock synchronization measurements. The Sigman 
calibrator continuously corrects for cable variations so that 
changes within an observation are automatically removed. 

In this appendix, the tone phase is specified in terms of the 
various effects that enter the calibrator signal in transit from 
the clock through the instrumentation and then through data 
reduction. Since the analysis is based on an "ideal" model 
for the system, extra effort might be required to account for 
deviations from ideal behavior. It is beyond the scope of 
this report to assess all complications. 

The model for the stopped tone phase is developed in the 
following sequence: (a) from the <t>-cal clock to the injection 

point, (b) from the injection point to the recorder and (c) 
through the data reduction procedure. The simplified block 
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diagram of the instrumentation shown in Fig. 6 will be useful in 
tracing the signal through the instrumentation. 

First, for step (a), the phase of the signal at the clock 
reference point (4>-cal clock) discussed in Section IV and 
illustrated in Figure 6 is given in terms of true time and 
nominal frequency by 


^C — U C ■*" T-c) 


(B.l) 


as indicated by Eqs . 4.3 and 4.5. We have chosen to represent 
the calibration phase in this way since (a) the clock error x c 
is explicitly shown, (b) true time t is common to both stations, 
and (c) tone frequencies are calculated at the correlator on the 
basis of the nominal frequency co c . The signal will experience a 
cable delay x u so that its phase entering the tone generator 
will be given by 


$c - U C + x c x u ) 


( B . 2 ) 


The "ideal" tone generator converts the zero crossings of the 
"x-MHz" input signal (5 MHz for the Rogers version, 20 MHz for 
the Sigman version) into rectangular pulses as indicated sche- 
matically in Fig. 8. Repetition rates slower than x MHz are 
obtained by "eliminating" pulses. The passed repetition rate 
will be denoted by u>p. As indicated below, the tone frequencies 
are the harmonics or the fundamental repetition rate tup. 

In practice, there will be deviations from the ideal waveforms 
shown in Fig. 8, deviations such as delays, phase shifts and an 
amplitude distortion arising in both the tone generator and 
connecting elements. For example, a zero crossing is not pre- 
cisely at the center of a rectangular pulse and the pulse shapes 
are not perfect rectangles. This report will not consider such 
deviations. (It should be pointed out that a constant deviation 
from rectangular ity will cause clock synchronization measurements 
to be biased by a constant error that is less than the pulse 
width but will cause no error in geophysical/astrometric 
measurements . ) 

The output of the tone generator can be decomposed into 
its harmonics, which will have frequencies given by noip (see 
Appendix F in [1]). The phase of the nt* 1 calibration harmonic 
(tone) at the injection point will be given by 


*n * “n (t + T c " V 


( B . 3 ) 
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The small group delays from the input of the tone generator to 
the injection point have not been explicitly represented since 
they (or any other group delay between the 4>-cal clock and 
injection point) can be included in t u . 

We are now prepared to model effects on the tone phase intro- 
duced in transit from the injection point to the recorder and 
then through the tone-stopping process. After injection, the 
tones, along with the natural source signal, pass through 
various filters and mixers until they reach baseband within a 
2-MHz bandwidth. As in the case of the natural source signal 
(Eq. 5.5), the effects of these instrumental stages can be 
represented by lumped variables so that the phase of the nth 
tone recorded at bit time t^k for station k becomes 


' t 'nk <t bk ) - “n (t pk + T ck ‘ T uk > " “hk s pk " *hk * ♦lk < "n ) 


(B . 4 ) 


k = i or j 


where tpk is given by Eq. 5.6. Again, tpk is the true time that 
corresponds to bit time t^k at station k. Note that the delay 
x a k does not appear, since the phase-calibration signal is not 
subjected to antenna delays. 

The recorded tones are processed at the same time the tapes are 
cross-correlated. The correlator uses a best estimate for the 
baseband frequency of each calibrator tone to separately "fringe- 
stop" each tone in a bitstream (see Appendix F in [1]). That 
is, in effect, the correlator will subtract 
the nth tone the phase 


<k * ( “n - “hk> (t bk _ 


which, based on Eq. 5.6, is also equal to 


^nk “ (w n “hk ) (t pk + T bk + T Ik 


from the phase of 
t s ) (B. 5) 


- t g ) (B . 6) 


where co^k is the best (nominal) estimate for the heterodyne 
frequency, tj^k is the bit time at station k, and t s is the 
reference time adopted by the correlator. For phase calibration 
to work properly, it is necessary to use the same values for 
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whk a n< 3 t s in both the tone model phase (Eq. B.5) and the cross- 
correlation model phase (Eq. 7.1). The stopped-tone phase at 
station k becomes 


. iT i x .T , T 

A W = * nk - t nk 

which, based on Eqs . B.4, B.5 and 5.6, is equal to 

4 +k < “n > = - “n <T bk ' T ck + T uk + T Ik - ts* 
- * Ik (» n ) - A$ hk 

where the "stopped" heterodyne phase is given by 

A *hk = u hk t pk + ^hk ” u hk (t bk “ t s ) 


( B . 7 ) 


(B . 8 ) 


(B . 9 ) 


This puts the stopped-tone phase in the form we want, expressed 
in terms of contributing effects due to instrumentation and 
correlator processing. This analysis has neglected some subtle 
but important problems that can seriously degrade the stopped 
tone phase. These problems arise from unwanted intermodulation 
of frequency components in the calibrated signal with frequency 
components in the stopping function, where the resultant, pos- 
sibly aliased, frequency of the intermodulation product happens 
to be close to the frequency of the stopped tone. Such problems 
have been greatly reduced in the BLKII processor by quantizing 
the counter-rotation tone sinusoids to 128 amplitude levels. 

Another fact ignored to this point is that the measured tone 
phase, like the interferometer phase, has integer-cycle ambigu- 
ities. Ambiguity resolution for tones falls into two categories: 
relative and absolute. By relative ambiguity resolution we mean 
the removal of ambiguities in a relative sense between tones 
within a channel passband. If relative ambiguities are not 
removed, the variation of phase with frequency within a passband 
cannot be determined and correct interpolations in frequency 
cannot be made. To remove relative ambiguities, one must have 
an a priori estimate of the slope of phase vs. frequency, which 
means one must have estimates for terms multiplying u n in 
Eq. B.8. Those terms are the correlator reference time t s , the 
calibrator cable delay (x u ), the instrumental delay from the 
tone injection point to the recorder (xj), and the epoch offset 
( T b " T c) between the recorder clock and the 4>-cal clock. The 
maximum allowed a priori error on these terms depends on tone 
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spacing. For example, if the closest tones are separated by 1 
MHz, then the collective a priori error (3a) on these terms 
must be less than 0.5 microseconds. 

After relative ambiguities have been removed, there will remain 
in the tone phases for each passband an overall integer-cycle 
ambiguity that is the same for all frequency points in the 
passband. That is, the tone phase will not be corrected for 
"absolute" integer cycles due to instrumentation. Accurate 
measurement and application of only the fractional part of the 
overall tone phase, however, adequately prepares the fringes for 
coherent addition over time and RF frequency and for an 
unambiguous estimate of BWS delay. 


B-6 



APPENDIX C 


REFORMULATION OF THE FRINGE FITTING EQUATIONS 


Insight into the least-squares fit to the fringes produced by 
Tensor can be enhanced by reformulating the output fringes under 
some simple assumptions. The reformulation leads, for example, 
to formulas for computing reference time t z , reference frequency 
(o z , and effective frequencies for charged-particle calibration. 

First, rewrite the expression for the expectation value for the 
observed fringes in Eq. 13.7 in the form 


V i4> Rk <t n ) 

< R(x f ,to f ) >= S c Z~t c fcn e (C.l) 

k f n 


where t n is the nt* 1 value for Greenwich time tg and where the 
fringe phase for bin k at time t n is given by 


*Rk <t n > = “k At R<V + ♦b + ♦e ‘ <u k - “ z ,t £ ‘ (t n " t z ) “f <c ' 2) 


and where Cfc n , which is theoretically given by Eq. 11.2, is the 
true relative channel amplitude, scaled in the ideal case to set 
the overall amplitude equal to the correlated flux S c . The 
subscript n has been added to c^ in order to explicitly denote 
the time dependence contained in bu(tg). It has been 
assumed here that the correlated flux is constant across the 
spanned bandwidth. 

The model fringes in Eq. 14.1 can be cast in the form 


Z m ^Rk'V 

c kn e < c - 3 > 

k , n 


where the model phase is given by 


+ Rk (t n> = <“k • “z> <T s - t £ ) + (t n “ t 2 ) <“k^s “ “f> + *s 


(C.4) 
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in which A s , x s , x<j> s , and 4> s are the adjustable parameters. The 

form of the model amplitude c^depends on the approach used to 

normalize the fringes. The derivation of this appendix is valid 
for both approaches to normalization discussed in Section XI. 

To determine how the fringe phase is implicitly fit with the 
model phase during estimation, consider a fit to the average 
fringes in which the sum 

L f = X | <R(T f ,« f )> - R m (T f ,w f ) | 2 (C. 5) 

T f ,a) f 


is minimized as a function of the four adjustable parameters. 
Completing the square, one obtains 

Lf = £ [l < R >! 2 + l R J 2 - < < R > R m + c.c.)] (C . 6) 

T £• / 


which can be rewritten as 


- r e e f. 

x^,co^ k,k' n,n' L 


S 2 c c e l(<<>Rk 4> ' Rk ‘ > 

b c c kn c k ' n ' e 


+ A 2 c m c m e l( ^ Rk ^ k ' } 

+ s c kn c k , n' e 

" (S c A s c kn c k ' n ' e + c - c ' ) 


(C . 7 ) 


The sum over xf, cof can now be carried out by interchanging the 
sum order, which results in sums of the form 


S = 

CD 


l(t n t n , )w f 


( C . 8 ) 


co. 
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and 


S 


T 


i(w k , -w k )x 
T f 


(C.9) 


As mentioned in the text, the least-squares fit to the fringes 
will be applied to (xf, cof) points in the vicinity of the peak 
point. For an accurate model (e.g., accurate values for the 

relative amplitudes c kn and fringes coherent over time and fre- 
quency as assumed in the model fringes in Eq. C.4), one would 
expect the parameter values estimated from a restricted set of 
points that include the bulk of the power in the peak area to be 
the same as the parameter values estimated from all points, on 
the average. (One can show that the partial derivatives of the 
fringes with respect to the estimated parameters reach maximum 
value in the vicinity of the peak so that the peak area dominates 
the solution.) Thus, for the purposes of this analysis, the fit 
region can be extended to cover all xf - ojf points produced by 
Tensor. Given this complete range, the DFT phaser sums in 
Eqs. C . 8 and C.9 are proportional to Kronecker delta functions: 


N -1 
u 

S = V e 
(0 / > 

k=0 


2-rri (n-n 1 )k/N 


u « 6 , 

nn 


and 


N -1 

2-rri (k-k ' )n/N 

S x = 2^ 6 “ 6 kk- 

n=0 


(C.10) 


(C. 11) 


The double sums over k and over n in Eq. C.7 each reduce to a 
single sum, which gives 

Lf “ S °^ n [ s ^ + a * ■ 2ScAs ° os A ^ (tn> ] (c,i2> 


where it has been assumed that the model amplitude is very 
accurate (c kn = c kn ) and where the residual phase is given by 

A» ok < t n ) = WV - Yc'-VV <c 
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in which the true phase is given by 


*c < “k' t n > = VW + VVV + VvV (cu) 


and the model phase by 


ur m (co, ,t ) = (u. — co )t + (t — t )u), t. + <J> (C . 15) 

T c v k n 7 v k z 7 s v n z 7 k 4>s T s 


If the model parameters are close to the final values and if the 
instrumentation is coherent in both time and frequency, the 
difference in the true phase and model phase will be a small 
fraction of a cycle for all values of k and t. Expanding the 
cosine under these assumptions, one obtains 
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F 


a 



2 




+ S A 
c s 




<V 


(C.16) 


where terms higher than second order in differences have been 
discarded. (It should be noted that the phase nonlinearity over 
frequency introduced by the ionosphere will cause a small amount 
of incoherency across the spanned bandwidth. However, for 
current spanned bandwidths at S band (< 85 MHz at 2300 MHz) 
and at X band (< 400 MHz at 8400 MHz), the deviation from a 
linear phase model is sufficiently small to make this cosine 
expansion a good approximation, even for the worst-case iono- 
sphere. For lower observing frequencies and/or larger fractional 
spanned bandwidths, the approximation should be reconsidered.) 
Thus, under the assumptions stated above, the least-squares fit 
of Eq. C.16 is, in effect, equivalent to two separate least- 
squares fits: one to the overall fringe amplitude (the first 
term in Eq. C.16) and one to the implicit phase values for^each 
time and frequency point, with relative weights equal to c^ n 
(the second term in Eq. C.16). ^ 

Given the least-squares form in Eq. C.16, one can readily show 
that the error in the output phase is decreased and that the 
correlation # between the estimates for phase 4> s and for phase- 
delay rate *<j> s is reduced if the reference time is given by 



k,n 
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It has been assumed that ( the coefficients Ck n are nearly constant 
or factor (i.e., c^n = c^ c^) • Similarly, the error in the 
output phase is decreased and the correlation between the 
estimates for BWS delay x s and phase <t> s is reduced if the 
reference frequency is given by 


co 


z 



( C . 18 ) 


(After being calculated according to the last two formulas, 
these two reference values are imposed on the data through the 
counter-rotation in Eq. 13.6 and in the formation of the fringe 
fitting model in Eq. 14.1.) Given these definitions for t 2 and 
co z and the above assumption concerning C]j n , one can easily show 
that the least-squares fit to phase in Eq. C.16 is equivalent to 
the following expressions for the estimated parameters: 
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<(co - g) z )y c > c 

«w - C0 z ) 2 > c 


(C . 19 ) 


; <(t - t )u>y > 

T^g = Z_ — <L_c_ (c . 20 ) 

«t - t z )V > 0 


= <»o>c 


(C . 21 ) 


vhere the averaging brackets are defined by 


<f (w,t)Y e > 0 = 


C kn £ ^ w k ,t n^ V'VV 


E 

k,n 


(C . 22) 


'kn 


and where ^/ c is given by Eq. C.14. 

This analysis has provided formulas for computing the reference 
frequency and reference time used in data reduction and provided 
a method for computing how a given phase effect maps into the 
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extracted observables. The validity of these equations, however, 
is dependent on the assumption that the model fringes in Eq. C.3 
accurately represent the actual fringes. For that assumption to 
be valid, relative bin amplitudes must accurately model the 
actual relative bin amplitudes and the observed fringes must be 
sufficiently coherent over time and frequency. As indicated in 
the text, measures should be taken to insure the accuracy of the 
model . 
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APPENDIX D 


DATA STRUCTURE IN THE TENSOR SOFTWARE 


The data structure discussed in this appendix is specific to the 
current Tensor software. The discussion is limited to a high- 
level description, since it is outside the scope of this report 
to treat details of implementation and operation. 

Fringes supplied to the Tensor software by the correlator are 
uniformly sampled in time but are sampled in RF frequency only 
within each 2-MHz channel. Sample spacing in time will be 
adjusted by the user according to the maximum residual fringe 
rate, while the frequency spacing will normally be 0 . 5 MHz 
between channel bins. For example, input points to the FFT 
along the frequency axis would number 800 for a 400-MHz spanned 
bandwidth, including the zero points used to fill the gaps 
between channels. However, the total number of points along the 
frequency axis will be 2048 after "rounding" 800 to the next 
highest value of 2**N (1024) and appending an array of 1024 zero 
points in order to double the sample rate in the output (BWS) 
delay variable. For this example, sample spacing along the 

output delay axis is (0.5 MHz • 2048 ) -1 ~ 1 nsec. 

Figure 5a illustrates how the four bins from each of the 2-MHz 
channels are inserted in the 2048-point FFT buffer for a given 
integration interval of length T a . These points are trans- 
formed by an FFT over frequency to the BWS-delay domain. After 
completion of the FFT for each time point, the output is saved 
for the subsequent FFT over time. Since the a priori delay 
uncertainty is usually much smaller than the span of delays 
covered in the FFT output, only a window of output points around 
the a priori delay value (at zero residual delay) is stored in 
the array memory, as illustrated in Figure 5a. For an a priori 
delay uncertainty of +250 nsec (+1 lag), only 512 of the 2048 
output delay points are retained in this example. Note that for 
a standard FFT, the output points corresponding to negative 
delays are contained in the upper half of the FFT output and 
must be shifted to their proper position. 

When all time points have been filled in, an FFT over time is 
performed on the fringes for each BWS delay. In analogy with 
the transform over frequency, the upper half of this array (in 
the time direction) can be filled with zeros in order to double 
the sample rate along the output fringe-frequency axis. The 
fringes produced by Tensor are a relatively complicated function 
of the delay and fringe frequency with a main lobe centered on 
the actual values for the delay and fringe frequency, as illus- 
trated in Fig. 7. In the case of stable instrumentation and a 
low residual fringe rate, the functional form along the fringe 
frequency axis will be approximately equal to a sinx/x curve 
that peaks at the actual (residual) fringe frequency. Along the 
delay axis, the "delay resolution function" will peak at the 
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actual (residual) delay if the channels have been made coherent 
by phase calibration. The sidelobe structure of the delay 
resolution function will depend on channel placement within the 
spanned bandwidth, the relative amplitudes of the channel 
fringes, and the coherence between channels. 
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APPENDIX E 


ESTIMATION OF EFFECTIVE FREQUENCIES FOR 
CHARGED-PARTICLE CALIBRATION 


The observables in Eqs . 14.6, 14.7, and 14.8 contain effective 
frequencies that must be specified in terms of channel-bin 
frequencies. This appendix derives approximate expressions for 
those effective frequencies and explains how they are used for 
charged-particle calibration. 

The effect of charged particles on the fringe phase is given to 
good approximation at higher frequencies by 


4> e (w) 
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where co is the frequency, k e is a constant, and AN@ is the dif- 
ference in the integrated electron content along the two ray 
paths of the interferometer, in the order j minus i. The group 
delay corresponding to this phase shift is given by 
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while the corresponding phase delay is given by the same 
expression but with opposite sign. The time derivative of the 
phase delay is also proportional to co~ 2 . This co -2 dependence 
can be utilized to essentially eliminate the ionosphere as an 
error source. If measurements of the group delay, phase delay, 
or phase-delay rate are made at two widely separated RF frequen- 
cies (co^, 00^) , then the observed values from the two RF bands can 
he combined as follows to eliminate the ionosphere delay (or 
delay rate): 


°ab c a°a + c b°b 


where O a and are the observed values and where 
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If all other terms in the observables are non-dispersive (i.e., 
independent of frequency), those terms will "pass through" this 
calibration equation unchanged. 

Since interferometry observables in a given band are often 
extracted from a wide bandwidth of frequencies (e.g. , 400 MHz at 
8400 MHz), it is not always obvious which frequency value (i.e., 
effective frequency) should be used for that band in Eqs . E.4 
and E.5. This appendix derives formulas for estimating effec- 
tive frequencies for the BWS delay, phase-delay rate and phase 
delay when these observables are extracted from a wide frequency 
band by the method outlined in Section XIV. 

The importance of accurately estimating effective frequencies 
can be demonstrated by calculating the error in the calibrated 
delay that results from errors in those quantities. The delay 
values observed in two frequency bands are combined according to 
Eq. E.3 to obtain a corrected delay. If the errors in the esti- 
mates for the effective frequencies for the two bands are inde- 
pendent, one can readily show that the associated error in the 
calibrated delay is given by 
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where x e is the charged-particle delay at frequency cob and where 
Aco a , Acob are the errors in the effective frequencies. This 
expression can be used to derive a worst-case bound on the al- 
lowed errors in the effective frequencies. Suppose a) the max- 
imum allowed delay error due to this effect must be less than 
0.1 cm, b) the observing frequencies are S-band (2.3 GHz) and 
X-band (8.4 GHz), c) the 0 . I-cm error budget is equally divided 
between the two bands, and d) the worst-case charged-particle 
delay is 200 cm at X-band (i.e., the maximum ionosphere at 0 deg 
elevation over the last few years). Under these assumptions, 
the maximum allowed errors in the effective frequency are 0.35 
MHz at S-band and 1.4 MHz at X-band. Since the spanned band- 
width can currently be as large as 85 MHz at S-band and 400 MHz 
at X-band, these relatively tight constraints indicate that the 
effective frequencies should be determined with fairly good 
accuracy relative to the spanned bandwidth. 

Approximate expressions for effective frequencies to be used in 
dual-band calibration can be derived as follows for the three 
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types of observables. Eqs . C.19, C.20, and C.21 in Appendix C 
can be used to determine how a given phase effect maps into each 
of the three estimated parameters. For example, by substituting 
the expression in Eq. E.l in place of y in Eq. C.19, one can 
determine the average change in estimated BWS delay due to 
charged particles: 
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where a minus sign has been included to correspond with the 
convention in Eq. E.2. 

To reduce this expression to a useful form, one must separately 
consider each normalization option discussed in Section XII. 

When no normalization is applied, Eq. E.7 is difficult to eval- 
uate in the general case due to the possibility of time variation 
of b n in the coefficients c$ n (see Eq. 11.2). Thus, the simpli- 
fying assumption will be made that all time intervals for a given 
frequency channel have the same bit count, which leads to rela- 
tive amplitudes that are independent of time (i.e., of n) . The 
relative amplitudes are still given by Eq. 11.3 but with the 
understanding that in Eq. 11.3 might change from channel to 
channel but is the same for all time points for a given channel. 
The charged-particle delay in Eq. E.7 can now be rewritten in the 
form 
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where the effective frequency is defined by 
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in which c^ is given by Eq. 11.2 with the appropriate value for 
b n from Eq. 11.3 or 11.4. The quantity AN e denotes an 
unweighted time average across the observation. 
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In practice, actual bit counts will not necessarily conform to 
the aforementioned "stable-channel” assumption. Given a 
reasonably well-functioning system, however, this expression can 
still be used to obtain an accurate estimate of effective fre- 
quency if Nt in Eq. 11.3 is replaced by N^, the average of 
the bit counts for channel k across the observation. 

In an analogous fashion, one can use Eqs . C.20 and C.21 to 
compute the average ionosphere effect in delay rate: 
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and in phase: 
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where the effective frequencies are given by 
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Since phase-delay is obtained by dividing the phase by the 
reference frequency (see Eq. 14.8), the corresponding effective 
frequency for phase-delay becomes 
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where io z is the reference frequency defined by Eq. C.18. In 
Eq. E.10 the "averaqe rate of change" of electron content is 
defined by 


T Z (t n - V 4 W 
AN. = — 5 (E . 15 ) 
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One can readily show that this weighted average is equal to AN e 
if the electron content increases linearly over time. 

As an example calculation of effective frequencies, let the 
channel frequencies be given by 

s-band = 2000 + (0, 45, 60, 75, 80, 85) MHz 

x-band = 8000 + (0, 10, 40, 120, 210, 290, 340, 360) MHz 

Note that, for simplicity, bin frequencies within each channel 
are not specified. This simplification will not affect the 
substance of this example. Given these channel frequencies and 
equal channel amplitudes Ck n , Table E-l presents the effective 
frequencies predicted by Eqs . E.9, E.12 and E.14. 

An unweighted average of the channel frequencies has been 
included in the table for comparison. Note that there is a 
16— MHz spread in the S-band values even though the spanned band- 
width is only 85 MHz. A comparison of the spreads in these 
values with the worst-case accuracy requirements computed above 
(i.e., 0.35 MHz for S-band and 1.4 MHz for X-band) suggests that 
care must be taken in the computation of effective frequencies, 
particularly at S-band. 


Table E-l. An Example of Effective Frequencies 


Type 

S-band 

X-band 

BWS delay 

2041.76 

8174.96 

Phase delay 

2057.29 

8170.07 

Delay rate 

2057.70 

8172.43 

Average (co z ) 

2057.50 

8171.25 
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The results of this appendix depend on an assumption made in the 
derivation of Appendix C, namely that the model fringes accu- 
rately represent the actual fringes. If this assumption is not 
satisfied, the preceding expressions for calculating the effec- 
tive frequencies may not be sufficiently accurate. As indicated 
in the text, measures should be taken to insure the accuracy of 
the model . 
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